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ABSTRACT 


Bng@oOverabploOnseresearchsapproachals proposed for the 
optimal allocation.of input/output buffer memory among files 
SbOReaROOeVarlousyiLypes sof sauxiliary mass storage devices. 

The objectives are to minimize the input/output time and 
storage requirements of computer file systems. A mathematical 
modelwonpea file subsystem is constructed with a Giscrete non-— 
MinearsOnjecuive uncUPOn and discreve Tinear (constraints < 
Cperat Ons rescarcl weciniguece sch as Magrange muleipiaer 
method and non-linear programming are discussed, 

An extension of the problem from a Magnetic tape 
file system to a file system resident on direct access storage 
devices is introduced. Implementation problems are briefly 
mentioned. Illustrative numerical examples with sensitivity 


analysis are included in the appendices. 
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CHAPTER I 


INTRODUCTION 


in Presources sl locattonginevomputer Systems 


fhieruser Otmoperalions rescarchelechnigueseto Solve 
Droutensod *eallocavionsor laree corporave resources is 
relatively well established. [60] However, the application 
of these techniques to optimize resource allocation in 
COmpDUGEr Systems hasInote yetebeen widely meported.| Whe 
Pap vere custom snOue mainiy to the fact Ghat computer 
Syovcemo Navmertosliaves Nou been wetarce Porporave Pesource 
involving multi-million dollar capital or operating expendi- 
Lures wL hat Cradit lonalay occupy the attention-of operations 
Resewme nets. 

Liss ubwe bvOnwis Channeling. Wnowever. 2s=syscemns 
becomeimore® complex and expensive. “Consequently, the need 
for eamere ravionals and sysvemaulc approach Toe this problem 
ise becominepoincreasine ly eapparent. —= A Poremnos berequirement 
Mirren ittulbapol veavron- of these powerrul  vechniques sis 
Dnes corrects ddentvitieation of problem areas Aan computer 
Syoceiis. ve ouch Problems may exist iin’ tite file’ design, Jeb 
schedulsnee or hardware configuration or a System. 

Atterta iriver Literature survey in The Next’ section, 
LHL Smee sis Wikia concern A.tselr withthe Tdentila capiaonseand 
analysis of one such problem area, followed by a discussion 


eripogsibtertrestmentsvanduramitications of thesprobdem: 
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The appendices contain detailed numerical examples worked 
out for a small™subsystem to demonstrate the use of the 
proposed "mer nods "er *treavment - “Computer propram Jistines 
(in FORTRAN) are included for a complete documentation for 


the*purpose of implemenvation of the proposed project. 


Dee -eLreeravure "Survey 

The ™publvshed works referenced in this thesis fall 
into two broad categories, namely OR techniques and computer 
systems performance. 

There is not a great deal of published work on the 
application of “OR (constrained minimization/maximization) 
Techniques vo *the problem ‘otf *optamal allocation of computer 
systems resources.-Chu [ 9] discusses the optimal file allo- 
Gation in a mulbewple*computer system." He "creats the problem 
of optimal sharing of files in a computer network. Walker [32] 
and Waters [33] both present basically the same approach to 
opuimize “‘biocking Pactors on magnetic tape Tiles using the 
Lap cance melo paler tet hod. "= However, their treatments are 
limited to tape files used in an individual program. Woodrum 
[36] develops an analytic queueing model of floating buffers 
LOMmIinimiime: progceanm run time. “Hiss modelvis a semi-—Markovian 
process. He also makes a comparison between the performance 
Ofeimiaivirieepumreringsand=the Ustal double wburrering., silver, 
Curalsieor haves developed a Guantitative rule Tor automatic 
allocation of systems resources in a multiprogramming environ- 
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On the effects of component interactions on systems 
performance, Hellerman (19) analyzes the effect of overlap- 
ping processes on throughput in an idealized multiprogran- 
ming environment, and Chang and Wong (7) discuss channel 
interference.” File “organization methods ‘and itheirtperform- 
ancerware sarscusseu, Dy Abate Gteal {1), Chapin (8), Collméeyer 
(V1), and, Eowe (22-97 

Many other studies include various aspects of 
computer systems performance optimization, but are only 
marginally related to the use of OR techniques in constrained 
Minimization/maximization problems. They generally fall 
into the categories of queueing analysis and dynamic simula- 
ELON. 

In the field of OR techniques and algorithms, on 
the other hand, there is a wealth of published works on 
constrained minimization/maximization problems. Only a sub- 
set of the most relevant works are cited in the bibliography. 
Of immediate interest to this thesis are Rosen's gradient 
projection algorithm (82) and Everett's treatment of a gen- 
eralized Lagrange multiplier method (46). Other key refer- 
ences include Dorn's extensive survey of non-linear program- 
ming methods (72), textbook treatments of non-linear program- 
Mange yeAbaddcmO2) had lei (o/s) ahunziee tealia (729), 
Mangasarian (80), and Zangwill (87). Book treatments of 


integer programming include those by Hu (100), Muth and 


pean aie ine ge Sak. tenia 
siie¥ien Shedd hav ahollven attextangro: efi -exnorsYresmd Uy 
Seeeatind (2) degmHP’ C0) fo to edodh ys horeedeth orm sod 
«(28) swod bos , (£8) 
io. etssqau cucitavofulga) enibute 1adtto° ye =| +e as 
yino sta ted ,nobtesinitqe eonsoxolieg efeteye: tos 9gMOm 
ieeintrenss Gb soupindzs? AO tu seu edd of Lerater (i ianigzam 
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sdg2, A ind. .auolgoty dolgnrieixsa\aoigeieinia beniarteses, 
miqargolitid 601 ni boris sis eXttew taavale: recog of2 Yo tem 
toaihers 2"neeet sve ziaedd ata’ 09 gest814! ick ites 20° 
=“Pey.s To tramsaerd, <'T2e7url ne (28)  adsirzoals noir Qetowm 
syete:. 407 sant .C9N) bodsemwsilgtetum ogervgs! bz tage 7 
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“WaTZONG Wwonltenon Io stnomrnoxd Foodtevs ,( 15) vbodisa paim) 
pa, . AP dee taal (078) yotbel yice) oibeds od gite 
Ae Sensations Hook 0%) VKiwgdiad hie (28) erreesynad 
ths a Com a WS seit shun: antentnaery sauna 


a a: a — 


Thompson (105), Tonge (108), Wagner (61), and Weingartner 
(110). Survey papers on integer programming methods include 
those™ by Balinsk? (90) and Beale, (91). The problem of 
SUD=ODtImtzatiOns1s examined by Hitch (50) .) [he whole 
question of OR implementation is fully discussed in a book 
on this subject by Huysmans (52). 

dna see ce ta turemsuivey Con eCOMpUuLeY 5/5 COMSso} bla 
Migation nds Okemethodse i senecessalitvesketchy, andecthe 
reader is referred to a more comprehensive bibliography at 
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FRAME OF REFERENCE 


Cpl teorbanc Comm Auprodcch 
inehis paper jon thesdesien of the "THR'= 


multiprogramming system, Dijkstra [14] introduced the concept 
Cmenterarcti Cale testing Ort aucompilex system during its 
COncUrUCGULOnN. sHeaclaimmed stheacv wo syouvem 2s Soundiy constructed 
Dil lees OLOUr TT ynLes ved shOmero Me possible Cventualities vat 
Sachse raccny Olemncreasing —complexiuy..) [his hierarchical 
epproecn Wwiisiebe adopted here in establishing 4 frame of 
reLDerence 1n which system performance Gan be predicted with 


some measure of confidence. 
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C.e, ee nerdWwaree cont raurav. on 


In general, the hardware configuration of a typical 
econmercial davaliprocessin@anstallation feomprises a central 
Processing Unity.(CBU), several cores boxese of main’ memory, 
one or more selector channels, a multiplexor channel, several 
magnetic tape devices, disk devices, and other sundry 


pemlpheral mMevilces. 
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Main memory 


(Program A) 


Channel Channels 


binifer 


Higwecrl “hardware eonrisuration™ on "a typical’ coumercial 
Cave sUrocessineg cenver. 
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To a large extent, the input/output (1/0) control 
ynits.determine. thes characteristic manners. in) whichs the 
different, tbypes of sfiles, are accessed... The» unit-recordefiles 
(e.g. card reader and line printer) can only be accessed in 
Siperesuntts of, fixed sizecarecords a jThus-there dsqgnoal/0 
parameter pertaining to these files that. can be easily 
manipulated to improve the system performance. Furthermore, 
Paese Uni c-Krecordslalcsmand devicessarcad i tected tosEne 
MUlgipiexor channe! sawhicheisestructurad ly. andviunetionally 
unepnnecledstOuthe hiph—speed peripneraledevi cess such as 
Giskeandstapet au hateis*s Chere is nor mutualschanne) inter- 
LekencenOramemory COntventAoneas such. betweensthenhighn—speed 
Pudge sow-cpeeheDerapheraladevVices; ekanelichteoteathissfunctional 
Separation, one scan _reasonablysisolate and) remove Ghat part 
Oiee hes SsVStempand Lis Pequirements -associaved with, theslow-— 
Specdsaevicessl romannOverallesystemecensidcravion tier the 


BucpDese Oils eptimi zing while design and system performance. 


244 .tely HardwarenParameters 

Bach hardware device normally performs with a 
eharacteristiespeed and: load capacity... These. constitute the 
primary hardware parameters of the system. The ones which 
are Of Vinteresteineeneapresent Convext are the channel buffer 


capacity, SheaGecerdingedensitreshandsuransier rateswoL 


bythe channelabutfer capacity umay be, considered tobe, the 
longest chain of data that can be transferred by one I/0 
Somme ws Mees On theMmlLBMAS60 At 2s. 32 KB.) 
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taper anc disk tevices Une rotational ands seck times of disk 
devices; “and =the inver—-block cap lengethyand, svart—-stoo time 
of tape devices. Other secondary parameters such as I/0 
error receverVerates-willebe-considered,—butewLll not be 


incorporated into the mathematical model of the system. 
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iM Orderevo wOGuUS OUlsau Lenny Lom On LNG slLdchiat ica 
PuOn Ole es pLoOUleieanGmolserVe (the primary efnects wor she 
DrOpOsed UReatieny,, wel USmslitsteCOlstger The Ss impLes™ 
operating system -- a single-stream batch processing environ- 
Welw. | patel romeune MeCescary core—res dele tOnuroL 
PEOeramSs, Ciere as only Cone application program residing, in 
Une -CiLaresrematnitte Coker Memory.) Luts Sxecuved Miva Ov 
COmplLeu Lom Dehore, Ene nex» JOD Le Drought Imvo core memory. 
Wiauete.  vthiere !oeno OVerldp ot CPU and T/0 operations 
Bepweecll dil iereny programs. 

RIV taperr less wine ave fixedes. Ze records thas 
Gonnve op rocked.., Disk miles; Well iOveue, wconssoered In Chis 
Syctemunittally., since tney involve many more paramevers 


GQuUeHTOrthe ditierenv methods of organizavion and access. 


Beles SOOmMUWOL eRe alamouel s 

Vie oar eliiOMmesLoes Ol MaLt Memory (Can pe a lXed met 
system generation time, or dynamically variable according to 
users' requirements. The maximum number of alternating I/0 


buffers per file ie determined by the I/0 control program. 
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IDENTIFICATION OF PROBLEM 


Kah) | lipmey-Veyets mops Blocking on Tape File ysystems 


Im a typical datasprocessimoinstallation of a firm, 
the majority of regular production jobs involve magnetic tape 
files. These files, by the nature of their recording medium 
and access mechanism, lend themselves to the various sequen- 
tial access methods. 

It is an elementary fact in data processing that 
such files can be read or written in a shorter time” if the 
individual records are grouped together in blocks of physic- 
ally contiguous records, separated by inter-block gaps (IBG's). 
This physical blocking of records reduces the total number 
of IBG's on a’ file. Thus the larger the files are, the more 
substantial *wrll*be the tape “savime duesto {ile *compaction. 

Figure 3.1 shows how blocking affects the relative 
amount of tape wasted compared to the length of tape actually 
containing records.’ It is quite apparent that the bulk of 
the tape contains empty IBG's when the file is unblocked 
Gi. 6a DLoCckingstacton =.l )s saWben ab locked by a clactonsoLetwo 
CBF=02,) jathere tses tall considerable stape was tace.., However, 
ltmsNOmLdebDeenoteder LOMathcac ra pismilieag Cures) 5.2. and. 5).5 
thatetheschancewot  Bretrom Isto Zi yicldssthe greatest 
percentage reduction (almost 50%) in total tape length 
requirement and I/O time per file. 

2 The time referred to here is the total I/O time 


(including overhead) rather than the total channel transfer 
time, which is independent of the blocking factor used. 
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Direction of; tape notion 
Blockineslact ores. 


BF=2 logical record 


Hie, 3.1 Hitect of blocking —(S0-byterelocicad records) 
on IBM/2415 models 4-6 tape device, 9-track, 1600 bpi, 
no overlap. (Actual size) 


aks (@nmengmictewiIN AS Gicons™ ieb al stent \Claiteks. Rekha! be blocked, the 
CMpire ii Se Vshouldscompr se4Oneso Loch wall owine 1G GO pe 


processed with only one 1/0 command. “In oractice, however, 


S: 


tis Le enero ecoctD Le sores trap les Hence... some compro 


Mise (has tesbe weaohed “invdetermenane enc sopcimalsolocking 


hacveor subject Osvue Civen cConstrarnius. 


3 In practice, the size constraints of main memory and 
channel buffer capacity impose an upper bound well below the 
sige of most files, Furthermore } Statistical 1/0 error rates 
involving reread/rewrite of entire blocks become a significant 
Cone Lderaertil Onvas sche tlocksi7e becomes, too large. 
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BLOCKING FACTORS 


2 TAPE I/O TIME PER FILE AS FUNCTION OF B.F, 
(TBM/2415 MOD 4-6.TAPE DEVICE, 9-T, 1600 BPI, NO OVERLAP) 
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BLOCKING FACTORS 
TAPE LENGTH REQUIRED PER FILE AS FUNCTION OF B.F. 


(1BM/2415 MOD 4-6, TAPE DEVICE, 9-T, 1600 BPI, NO OVERLAP) 


FIG, 3,3 
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BLOCKING FACTORS 


FIG, 3,4 TAPE I/O TIME PER FILE AS FUNCTION OF B.F. 
(IBM/2415 MOD 4-6 TAPE DEVICE, 9-T, 1600 BPI, NO OVERLAP) 
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BLOCKING FACTORS 


FIG, 3.5 TAPE LENGTH REQUIRED PER FILE AS FUNCTION OF B.F. 
(IBM/2415 MOD 4-6 TAPE DEVICE, 9-T, 1600 BPI, NO OVERLAP) 
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3, lt) FETIPSecrs ony /OeTime and Tape Length 


In general, the inter-block start-stop time per 
block is independent of the block Anne Thus the time 


requ ned to read/write a block ws piven by 

treaeortecex-c (864) 
Hence, the read/write time per file is given by 

De= KExR eset Caxts) (362) 
Simibarhystthe topaletapethengsthyrequired perrrilesis teiven by 

a= K pee CGetD=xoB) (3x5) 


where Ge= B/Oetinerpertbicck 

Ep=nb/AC oime iver iiiie 

i =ulapemilenath (per tiishe 

Dae=aTape Length perscharacter 

Co=al/Octime per character 

S = Inter-block start-stop time 

G = Inter-block gap length 

Kea NO Otepockssper siavle 

B= No. of characters per block 
ul PatsSeLSher laste ess ump clon in our proolemerdenuiitess won. 
[im bach it Varies over a range which depends on the nanuiac— 
furer and the model of the tape device, as well as the file 
Sreanimab.on and the processing time of the program. yin some 


instances, the tape movement need not come to a complete stop 
in between blocks. Here, the mean start-stop time is used. 


ReomeeduetblCnses.2 and 3553, "the total 1/0 tame T 
and tape length L can be determined with respect to the 


Broeckastzens? 


Sv PErecuctl on LoOrely/Oulimevandmrovernead 
Promirigurevys.2@one lean readily sobserve’tharethe 
MOStPsligniftveanh variations imPLZOnGime. Mvecinathe rance of 
Di¢gewinesiactors#betweenS ll CandPvoOresivsis alsevevident eircom 
thepwrantiytoilctwvessinuiicurens.4othatm Lor wanysparticular 
dirirerencetin Dileckingpractors , Geachabsoluvetvariationyim 
EZOStimeiisedirect lyVoroport tonal te theaktile size jnwhereas 
Pieaperncentace varpalloniretindependens tot theeiiie suze: 
Tiemireductioneef ctheminumbertor EBG's"on, aetape fie 
VReLas Prinmakiihyearedquclioneot athe CoObalbeanter=bilock 
Svarv-svopmlime.) Concomi pantewsunuthis [shiunermreductiontorl 
i7© overneads whichtconsitsts mom. 1/0 oommands) and Ssupernvisory 
routines: for handling *i/0 interrupts*(and programmswitching/ 
swapping in multiprogramming/time-sharing systems). 
Themlabiereareesecondanry lonltentiary .6f tecrsuwith whichiwe 
SheiiineoumoePeOncernedhin -OounT basi eriranenor Gnenerence. 
These effects will besdiseussed more fully in aymore complex 


Prame ofereferencesat a later stage. 
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Svuet,. ereape Cost savings 


hate wet Eworcvhenotinemiromel! curems.Satbaar most 
transactions filesporelesssthan aboutel07000 recordsr (of B0= 
byte logical record length), when blocked, can easily be 
stored on mini-reels of 300-400 feet long, as compared to 
full reels of 1200-2400 feet long. When applied to the 
entire’ DP systemfor allarge corporarion which uses 15000 to 
py OUUmreels of Vapes thls can @represente a substantial: saving 
iether cost, Of Mbapewacouis 1 Tons) for an instal ation, 

Two basic approaches can be taken to reduce the 
Capewrequirements. »One can svtoresvall small fries! on Separate 
Mine=Ceeis,, Ot combinevseverai small Thles into one muici— 
VOTQinCw UIC Ol om Ui Peel emcee pS NOUN Cm DemvaAKel) GOLUSe ie 
batter lacproachriTceonjunecuroOn with the job scheduling, 

SO ica Tne wsmalie Tiles ioOrvwcenseculvive (obs tare Ss voredeon 
the same reel in order to minimize the operator's tape set-up 
time.) in a single—-streéam batch processing environment , 
usualy i ound Ane smaller installations, an excessiveyset—up 
tine, can become an important factor in dercrading the through— 
PUcwOtmeunecmoysvem, = Hence Ther mini—-reel approach mene pus 
such a system at a disadvantage. Ina larger multiprogramming 
environment, however, the set-up time can usually be imbedded 
iam chiesOrocess ine UiMmer ol LOvher ODS monrough inv errniod 
scheduling (such as HASP), by which the operator is notified 
of set-up requirements of impending jobs in advance on the 


Revyooard console, ~in such ae situavlom, the mine —reel aprroacy 


¥ a ope ian pes eats 
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WoUbasprey lide smorer tex ips apy eo escheduling without surferine 


The above-mentionedvdrawback. PNormaliy, a combination of 


these two basic approaches would be adopted in order to make 


full eusemeonetnescharactermsracs Of “each apo lication subsystem. 


7 D a te Ae ean _ 
en ee Je oe Se 
- P . - + an 
a 


oe one, 


7 


es) 


CHAPTER IV 


TREATMENT OF PROBLEM 


4.1 Methodology of Treatment 


Dee we cinciples: of Opuimization 


Optimization of tape’ I/O time and length require- 
ment can be @chteved through the mathematical determitiation 
Cuanpeochano Dae nOnsattisans the caleulus, Of variations and 
Hou -atieas Pe prorrammnime metnods., “iitprincinle, cre tries te 
minimize (or maximize) an objective function subject to 
Cores ino system Cconstrainyss 

For? esoeorram using two or mone tape fales. ca 
Certain upper bDeunteonkthe sizes or theyl/O0 burrers is) imposed 
Dy system constraints ~ (Ue.2. job pareition size, program size, 
@hannel butfér capacity ,suavisuncal 1/0 error rate, etc. ), 
dependent on what hardware configuration is used, (e.g. size 
of main memory, number of selector channels, etc.), and what 


software support is available, (e.g. 


© 


Single-stream batch, 
multiprogramming, etc.). 

Ths weper pound Limiipeepoene one nov duole 7 O 
Dui ter sizes. snd. oheireagerecate Size pine obj J ecliver as uo 
find an optimal arrangement of blocking factors which would 
Minimize thei total number ofecd BG" "son thermape Tadies subject 
to this buffer memory constraint. 

One might ask at this point whether the minimization 


Or tne number of [BG's is The imest Suitebic <riterion Gorchoose 
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for optimizing: the isystem performance’. © ihe ichoice sofian 
COPEC GIN ere cura CuuelLatiparcto mibhericonrecnWdentification 
Sno Ursavienerol the problems Accordingly, thls quesiticn 
will be dealt with more fully in section 4.2 on the 


FOrmulavlon. Of wLNe Ob ;ecei ver function. 
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To optimize the I/O time and tape length requirement, 
5 


we minimize the non-linear objective function 


ae . Fo 
eS eames + x Gee) 
4 n 
sub, e6cypetoutLhenlinégar memory cons traintvequation 
copy cad Ghee ae i (4.2) 
AG 2an)'s gei525 sf 
where Zh =, OvalenG  MoriplocksHCieGia) *pracessedper tiob 
M = Total core memory available for tape buffers 
= = Blocking factor for file j (records per block) 
B = Size of file j (records per file) 
a = No {fOr alternare wapecebuiters sper fier] 
Re = Record size for file j (characters per record) 
So = Slack variable 


In equations 4.1 and 4.2, the parameter a and M 
Gan be Ancrementally varied To observe their e€fiects on she 


SOUR LOns sor ole ce: Various techniques can be used 


5 This model .can bé reformulated to consist of a linear 
Objective function subjectvero a*non-linear*constraimg. “Buu, 
tiieeneraly the former model is” easier to solve with the 
available computer algorithms. 
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to find the desired) solutions. The: mathematical techniques 
applied) to. solvesthesmodelacanabes ereatiy simplafied by 
OUSeCrVing the Pact that the discrete objective function z 

is asimple jewe bi-pehnaved funection; which can be treated, as 
STSMOOTh, —CONCINUOUSH IT UNnCGULON With No Sincularitve wach none 
Hegouolerceri On westOr sinmplicuLy. or ~aporonch.. Lie Mepnod som 
Lagrange multiplier [46,57] can easily yield accurate results 


by hand, calculations or with the;use of a\short program. 


4.1.2.1 Lagrange Multiplier Method 


Taka neupeartiadecerivaravessofane eandystivain 
equations 4.1 and 4.2 with respect to X, and equating to 
ZerOuLores Wlad , WeLodbaline the ordinary rpartiaiecitiercnyial 


equations 


AO) 2 nee Gea: eet AME i ee tal (4.3) 


Vere Me Soule barranven muy ior Perr 


6 Meassumes nO DnYSs1calescipniticancewin. OUT cOluuton, sas 
We are only anverested in solving @ror ths Voie, the Dlocking 
facvors. However, an amportant economic interpretation is 
apvuacned to A, Which represents the incremental reduction 

Mai une oval number of blocks processed when tne vovsly aval l— 
able buffer memory is incremented. 


2 


= - = . 
rhe beat Sar’ 


a 
a Of cae 
4:2 vine, 08 i age 


‘Saepditape fav testes SYP. 3 A ie 
4 . — i 
ya . 194 ipa! Pouca eC 108 ha Ts 14 


: © ; 4 A 
nt Seite wa 9 -sifi sng Gade 47@3 eés. 5 THES 
Atrrnbt Apendoattaw ,etacte a me 
ta - ‘ _ 
“oe ; ’ | int?e _ ri , ay) * Aton 
é i ae | 8 t : 
e ' ’ ys] yy a OF ’ y 
’ 4 = de . at | 
i 
—— 
— _ pis een a 
4 LA A t) i j ‘- y Oe 
; pl " Dw ese «=? #8 
j oe} : ’ ia 5 .) ees ae 
s 
.) =) 
= 4 - = 2 
7 j 
tT] ' i 
[ i] ut 
‘ 
_ 7 
> iar a 
a ; : ‘ a - 7 
; 7 . — - - _ 
gaolsulus pl 4) an oa aFig ngd™ - aes _ ~~) Sri &e a 
- 2 » | =? i? 7 vie » SV im 
ef9 e013" 40° 5. OE. 30) Ot cere tay rr 


te » 0@ ve 


- 23 


From (4.3) we obtain a set of n homogeneous 


Simultaneous algebraic equations in ntl unknowns 
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La eh = = 

2 + ask, r Os i= our eee erie wal) (4a) 
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Combining (4.4) with the original memory constraint equation, 
Wemcane Vil Stisoly cca Weebratocal Vy OTaarele elma a! She pele atieic 
Sue Be and wrecondy size Ry are Knowmrconstantsiitor ach file,’ 
and as and ee Metered ecwsion variabies + 

hey should be reniembered thats tne objective function 
ana Variables in whe pivsitcal system ican only take on positive 
iivereravalues. | Tnewmmathnematicallyfoptained Optimal Solution 
Can thengve roundet=orr Neuristically and tesved for Teasipi— 
LiGyeanOmopu LNA ty COecivera physilcally optimal or mear 
Opolmalesolutitons FP The unused, port von sor the "available buiter 


Wemory grverra solution 18 rounGed™@ oF pis Simply taken up by 


the slack variable “5S. 


i ineeeneral, the file size ls vends vO vary Over a per lod 
GPavime = in some cases “considerably “Under those cir cum-— 
Stances, Bs can only be approximated by its mean size. 
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4.1.3 Global Optimization of Entire System 


When the application programs involving tape files 
ino datas processinge departments arexconsideredmasean interrated 
Syoleom, ere When someair iwles  aregintertmedarcdethrougchs more 
Chenwene Hrogram, the optvimizaciontored lem takes ontarnew 
aamenson of complexity. 

Here we shall define a ‘shared file' as one* which 
Serves ase i npuLriTOnmore. thane oner program, and a "transmitted 
fote'las oneswhich as thet outputs from. one: program’ and® the 


Lops LOmanothersprocram:. 
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Consiner TheyproplemmthisewaysurAstapes files} is 
sriaredsbeuweenrerorrams A and Be “if its bilockine factor as 
calculated subject to the memory constraint equation 
associ apedawithiprogram,A,#(hocalxoptimizgation) .42t is 
assigned a value as Similarly, using the memory constraint 
Squads sone lor proecranm bh. Aes is ODcatned. foince ws terse nehact 
thewoeme Tile being considered. it must necessarily be 
asslened akunique blockine®factor o- Coleutaved sSubvecu to 
Doth the, original two constraints, which become inequalities 
(less than or equal) in order to yield a unique solution. 

Ta as important to readize thavlin acomplex system, the 
optimal x, MEeCdANOEBneccessarmlynwe Kay or ane oretheirvaverage, 
andkanlarbitraryrassignmentmor * wouldwentatipereadjustment 
Opepnecsoris Derkettetncepapeminiles in the subsystem: Often, 
SUCH Geperturnbavblonson agcomplex@systemtis surrficientyto 
Hniduceha sieniiicanvedevtacionvoreasy stem i pernormancesrrom 

Bie opormaw. PTH seisswhyeopoimization slechniques@such@as 
Nen—-linearsprocramminghare necessaryaro ensurewaetruly 
Opuvmmalmsolutiont 

Let us consider the igeneral icase of a subsystem of 
meinterrelated application*programs involvingen tape fites’. 


THES torn OLeche: Obyectlve:, functloOnyZ remains the same vas Lor 


a singe fade cron. © Forecachsprocram, JtUperesiseaniinequalaty 


8 The objective function now applies ton the entire ysub= 
system as a whole, although, for the same number of-tape 
Grléeswoitpiteakes thevsame iiorm asginethetcasemoipe Locad 
optimization *or aasinglerprogram: 
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menoryeconstrainteteThusswerminimizeszesubject tovansetuof 


m memory constraints.” 
a5,R,%, +8 ct apa aearion cf Sy < M, 
(Ae 59 
ane le, "3 is Fons n i Sin n Min 
As as (a Oo aan 
where 45° oe Ky and Ss haves chersamey meaningwasvyansthesecase 


Of Tansingdesprocramveuhere M, denotes the maximum available 
Taper Dulieramenocryeperrprorcrams is 

whe solutes onmpos thnase problem cangibea round. by 
epumyingr certain non-linear programnminge techni ques  suchaas 
the Gradient Projection Method.[82, 83] These techniques 
have been summarized in a paper by Dorn [72], who surveys the 
qgevelopments anivthis fleldsuptrosthesearsy sixties: 
A Gheoreticalsexposition’ of theserymethodsywouddsbe beyond the 


Scopey Ops thist applicatton-orlentediaiseertation Sand nence 


9 Same Cpetnescoeritcients of AU. esLes ome! elenents ot 

the matrix (aR)s;5 will be zero. J The physical interpreta- 
TO nroOnouiie) Ls that notesgiletheviades witinbe used sagaay 
Gartlcular program, Thus in general, the larger a .subsystem 
is, the more sparse will be the memory constraints matrix. 
iiehiphveot ebhis, characverisv1lc. val eoritnms thataare sapetl a ley 
well-adapted to solving large sparse matrices (of the order 

of say 50 x 30) can be employed when dealing with large 
subsystems, (e.g. with 20 programs and 30 tape files). 
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Wi li not bewartempred.  sutiwoce 10 to say that given 2 program 
whieh implements ajsuitable alisorithm, a person trained in 
Ssucnetechnigques Can easily interpret the numerical resulve. 
However, to set up a mathematical model to represent this 
PoOviemswouLdecequire avskiitedMansiyst.of commurverrsysmens 
DOseValuate thesrelative sligniiicancetof numerouswinterre lated 


Der anegers 


4.2 Formulation of the Objective Function 

MRESob ective hunet ion zacam rebel conssderedgasiGne 
total number of tape blocks (IBG's) processed (read/write) 
Dee wllcs SG ressuls ys CeMeInmOne ri cycle tar ASsubsyecemeis 
defined here as a group of programs (and their associated 
files) in which egeupnogram tetreillatedivo atbeleast one 
other program ih the group through a- straddled file (i.e. 
shared or transmitted files collectively). In our present 
Conivexty aarunmeyclesoréatsubsystemmisidetimedtestune nun 
ey CLerorrunedieast efrequenvlywruntproeramrinithe eroup. 

Rreagetapetuisgeris straddled; say, -Overpprogcrans 
A and B, and A is run daily (5-day week) whereas B is run 
weekiy, WuneneA should? haveraGrelervivenwedentaone seston iin 
Tovalitns verthe anumben ofritape bDilockssprocessed in that 
period. Thus to reflect the relative weight (frequency) of 
the number of blocks processed by each program with respect 


to each file, the general equation 4.1 should be modified to 


- 
7 
7 : , 7 a * 
7 hae 1 ee ad 12 oar ph, “ae wiqears 


=r le} 3 gece 


a 
set eeh |, v reVia 


; “a 
vVisieromce Ba 


, . = 
» 7a. | = 7 
wy, 0 
71.1 2 < = ag 1 
2 2 ne } # 
é ' 7 Ls ’ 
miu « ant ; 
‘ ‘ 
' p “ ‘* ay 
7 + » 
ad 7 oe i] é if Ay 
‘ « i$ ma 
} * - € ‘ * 
‘ i , iil - 3 
) =i - bi ‘~ = 
2 ‘ i at Mia 
gs @ ie Vnete : i 
= P 
1 teri) sh  WeeeeeoT). cool o.8 yi ae ons, w 
: = fee 
3% ,@9a0c ner fib: ug? af’) Fol eu heats _ 2° eal} 


4 - _ 
tages eco me tage phen me = ineate 
. oY) 


ene te i nse i A Gphiwe Le on + thie 4 


- 28 


F F 
hy al n 
ZS ea es ae Se (4.6) 
1 n 
where i. Gh=b,25 s2nten) shisishetrelativesnumbercotatines 
PVlewiwees DROcesseasimthant tne tsubsystvem nun. cyclers tine 


HOcmIULaAs horn 4gcould Thenape Vsederosgesl culate the total tape 
I/O time for the entire subsystem over one cycle time. 
AMeOUDerEEIGrOr Which GContribuves To the number oF 
blocks processed per file is the number of reread/rewrite 
due to I/O error recovery attempts on the tape devices. This 
NeenOcle Geral tho hOriaranGacvOornlOnmincorporave iIntogtne 
Obj) ee vlyesrtunecth1 en swwecause OT el TSsmstaelristicatinatures 
Furthermore, the-sequence of motions which the’ tape goes 
through to retry an error block is device-dependent, (i.e. 
dependent on the device model and the manufacturer). In 
Srderetcadevermine tts <~fiecse , sonemwouldgneecdyccorknow the 
IAORSCrORnVe Lab Sbics TpDergyoluMesoOVvemnasperiod of time, the 
user/system specified error threshold, and the operational 
enavracrer hs tl Copgortthermtaposdeyvicecsemreor Ordinary biock 
Bigeee However ssphisaiacvor canmbesconsidered  nesincipile 


(for well kept tapes and read/write heads). 
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Pheiquestion has@bpeen®posed in =section™4.2@itas to 
whether the total number of IBG's processed by a subsystem 
Ls indeed thesmostesuitanlescriterion forlest imeating) incre— 
mental changes in system throughput and cost-performance ratio. 
Togeniswer thus sques tion serrectively ~ ones has tomfiirsricitall 


deal with these two notions separately. 


UM mir li MECostaCeiveria 

The meanime %or cost is eby novmeans self-evident. 
Lewcould Lrererstosthe dollfars=and-cents coer -Ofrethertotall 
NaenWetl we wane acouisittonseoL an installations, Tuecouldebpe 
ine verms of the amount -of -coresmemory thay cansabe removed 
from one section (partition, program) and expended for tape 
buitering Mireanother ly 7Orsagain ivecould mean-the quantum 
Jump In VOperavine personnelecost 1anvaddingtanother work=shirt 
TperneLCOMputer sCenver. P@some OL SGReSeecost Lune vvons meehu 
indeed involve mutually opposing trade-offs: 

it the Cost Ol Peape -acqUu Us ft rons srpeinge voonsidered:, 
DOCH eDhnesOpbiMivat lon “orethe tora lenumb ereonelbG Vssalse 


serves to Minimize the total tape Length requirement: 


Vicelyeod elnnoushpulyCriverda 


The conventional meaning of throughput is more 
easily defined, namely the amount of useful work (the number 


oftjobs) that “is “processed by the total -computer system in a 


ry 
~~ 
a 
<a 
é 


mew 


is: Ss pdaaie cedoe = aed tot 


al ast eae fae ita" 33 nigies tobe on? ventink i 
ga) nepeye ef “aaaiads teaczoem 


¥ r \ ° _ 
e a" a.) ie 
R C win ‘(joy l soethe nas rteaie she? Tense oe 
Pind oA <4) a @ ae are! Hw és * 
‘ = 7 ae 7 
é | - 
mh; 
.7 ‘ 
Jtny See Ss 
mt : ; 
a 
i Lory i ; ‘ } Pass ; nee as 


§ & ATS 4 ' = " ? 1 s La) ih 4 Sd 31 


7 a eo = s 
ce. set ay Senge 


mS 
a ea ss 
> 


— oo 


given amount’ of time. “This*derinition-is admittedly imprecise 
on several points, but nonetheless self-consistent as a 
Gencepcuel definitdon. 

UR SOUcGTOriMary Ob ;eculve ss 7a Nien tnrougnpurerave. 
then the amount of positive correlation between our original 
Opjective fPuncvion and une® sys ren eenroughnpuc tame woulda depend 
largely on our frame of reference, namely the complexity of 
the hardware/software configuration and application program. 
HnsourY baste” frame ol reverence, * the correlation as expected 
voy be-high’y/sFor® instances if a"207-réeduction Sf total™@tape 
T/O time is estimated from the optimization in a single-stream 
baven processing system which 13s "/0% 1/O8bound; the actuad 
reduction of throughput time can amount up to-14% (i.e. 20% 

OF 702% 

When this system acquires’ ‘more high-speed auxiliary 
memory devices (disks, drums) on the same selector channel, 
ehannwel= contention can become a-signiticant’ factor* to be 
HeCKONeCOPWLth = lis thet system crowsesanve; Semulviprocramming 
environment), the overlapping of *CeUa@ands1, © operationsrwils 
tAtroduce Turiner* secondary e€rrectsy on channel centenvi0n. 
Puriher Sop. svi catnons OLveune=@s Vox cll ecuch sace Uime—-smardvin, 
Wild™ introduce tertiary effects of overhead Trom- job *swapping, 
tines siacine memory paging, etc... ethe=inveracvion of ald 
these adcattiondl factors Will] rule out a sUraisnttorwara 
espimabion Of System performance On Tae Dbasiseor eheoproposed 


Opuimiuarion modeis “in such Cases yin order towbe more 
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determinate, sit willsserve well to employ a simulation model. 
HOWeVer, tUnat does not Ciminish the value of thie sootimizabion 


approach. 


4,3 Formulation of the Memory Constraints Matrix 

LHe MenorysCOusGraimmus Marri x 1s POormulatedendn 
Piece eoh sSsystbem restrictions “and memory avallaba lity 9) fen 
a rare case an individual block size (RX 5) obtained in the 
solution exceeds the channel buffer capacity, (e.g. 32 KB for 
the IBM/360), the associated channel capacity constraint 
(R,X, Ss max. channel buffer capacity) would have to be intro- 


ah 
aqueeds ‘he solution 2s then réecomputed So that “the max {X,} 


PELesewithin the upper bound imposed by the channel pGaree 
capacity. 

AMMOrPeSSyStemacvic approaches -oOf Course Goninciude 
all the extra constraints (RX, Ssitax Channel oudterecapacluy, 
ice en). HOWeVer,. thesinvrodyuction), of considerably 
more COnNStraints creates an unnecessarily unwieldy mode aes 

To evaluave the aneremental sdii ferences amone 
Variousralternabives in the optimization of blocking factors, 
One (Catia Vvary the decision parameters aaj and Ms This infor— 


MabplLon camuserve as an additional analytic Tool in estabii— 


SHIune Memory DarLioLoOne sizes) and | Obmechedules. 


10 ehor inetvance, in.a-subsystem,of l10yprograms sand 20 .tape 
files.,+the original matrix of 30 constraints would be expanded 
Om NOMIC OT St a Li iiss 
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Ao. Liebereats OhiManrpulabine as, 


OnesCanvespect. a Certain incremental gain in the 
Throughput of overlapped operations of CPU and I/O by increas-— 
Iie the number of alternate £70 buffers per tape fies” 
This = increase in a is normally effected at the expense of 
reducing the individual Durier sizes in a Constrained memory 
Saouavion. Ihus) this =gain, Gan be offser Dy tne Loss resulting 
PCOMeGhe Itcreasem= in. 1/0 Limemand Overhoad. This trade-off 
should be carefully weighed, especially when one is dealing 


with small blocking factors, whose incremental effects on 


I/O time and tape length requirement are much more pronounced. 


as Cae CCUS OF Vary ine M, 


Drees inerenental errectscolavaryine Une Tota ia Vall 
able amount of buffer memory can be used to estimate the cost- 
effectiveness of expanding a memory partition, or of reschedu- 
iiVipee  cerlLain jODeInlOra  larcers partitions lO. CeCcuce mit serun 
PeiCw tls paranelver se Ol maior Import niches ormilacion 
Of the mathematical model, since many aspects of sysvem 


performance differentials can be related to variations in Mo 


EE az; is subject to a maximum number that can be handled 
byetne [/Orcontrol system sottware. Ply sie mnormally ssslened 
auntie sot iheorae. 
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4.4 Important Considerations in System Analysis 


TOmde Si enWaneopuimalesySstem thate wil lenetalngite 
Revel Sorepertormance inwthettace of = futuresecrowtn and modifi-— 
Catlons to the hardware, software or job scheduling, the 
analys |) MUsv Necessarily mMmeakescompromises in Setting up a 
Mathematruca Wi mode i yroOsrepresenUethewGed ies duet tone euhis 
PGSSUDDOSe SS Uhawlinhormationgaboue such plans istavallabples 
HOWeVeGiewetis TRS maysclem redesien can be done in an economical 
Wayne cee eanolysuemay decide sho, make less val lowance {or 'ruvure 
changes set Toasts VOmsSvaie OP =DewLet COst-ellecuLVeneas 
ImMed Bateky.) OUCH Raumati syvsvem weuldsbe subject: ToOgmore 
heeqvenuredes.ems:. 

Although The primary coal en this project 1s) 10 
Optimvze: tape ficessystems, a system analysis would be incom— 
Pueve anc suMcound If tape operat lonst areca ccruts nized min 
PoOuaulen withouG considering how Chey Aanberactl Withy ovher 
DervomOre one ISVS vem. § SUCH GUCS U1 OMS footie =Cerrecesot alr wand 
t7 Oy overlap, Channel overlap. vand channel contention are 
Squally Aas importany as Uae mavnemavuca model, sal vnougie they 
may have to be solved by a stochastic simulation model without 
the power and elegance of a mathematical solution. 

By itself, the mathematical model (memory constraints 
matkixywonly describes the contention for main memory among 
Gielvaridousstape fides used inta program owhence P aus must pe 
lige, (eon Wine Ga Chuewith the eruherirelaved facvors Which are 


Summarized below: 
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Heey la kosenia al, Information 


head ol lowing cabercories: of information, are, consid— 


eredecccenl iA ampli (Om, AesOUund US ySLem analy sas. 


4.4.1.1 Hardware Configuration 
INMChey Ssysremearchitecturc. either analyvets should 
ascerlain whether the, CPUs s., accession the main memory is via 


a dedicated channel 


170 BUS 


MAIN 
MEMORY 


Or, Via. a .common..AO_ bus.. PERIPHERAL) DEVICES 


I/O BUS 


The former arrangement is by far more common and involves 
less channel interference. .in whe latter, case, which, is nov 
uncommon in small systems, memory access by the CPU may be 
handicapped by the malfunctionsof any 1/0,undt on tne sasyn- 
ehronous, 1/@ bus, om by an. inordinate amount,of channel 
contention among the peripheral units. 


Din addition, the number of peripheral devices 
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attached to veach*selectorachannel should be taken into consid-— 
eravlon . alpnougnV ite smditrecudcwloequanni tyerniseracrorean 
our model. <A system with two selector channels should have 
Inyecenervaltatlowerrrate lofichannelPinterferences _tASreneral 
intuitive assumption is thatusince large tape biocks monopolize 
ppeschonne Slorelongser bubereweribursts#of Otrancmisstoneicney 


benOPrtorcreatleupreaterpchannel contention, 


amet el oOftware support 


Pues cVeewO Opera cing ssyelemnUscOmnasmar GLroar 
HeCarlne “One cOresmemorys al locauion considerations, such as 
particion sizes WMrixedPorrvariabile,) Band) theirebasic unitmor 
boundary increments. The amount of main memory available 
PFOLeUSeCre DROP aliseIMy ar Dartdtienvas well as’ thel programisiazes 
(excluding the I/O buffers) must be determined. If a program 
fe wdavidedmintoeovertayincesezments Wehe avait@abteg burger 


mMemoryvemayabe Mimited Dyethe Size ofPrhbellearzcestasesmentr. 


4.4.1.3 Master System Flowchart 


A flowchart (or an equivalent description) of each 
application subsystem Showingethe aneerrelartons#o einddvidual 
procrame runs and tbheiryassotiatedi i leseasPingtepensable. 
ingaddition to, provyidingsanveéiiectives means, tovhiehniigng. the 
existence of bottlenecks in a system, its modular represent- 
ation also serves as an aid tovanalyzing different alternatives 


Stejobe scheduling. 
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4.4.1.4 System Audit Data 


The stavistical breakdown of audit data on system 
Usage indicates the amount of usage of “the CPU Sand 1/0 devices 
Ber eprogram=and per application. ~Percentace utilization of 
poese resourees and the) throvenpur rate can be calculatedy trom 
these cumulative data. The amount of I/O overhead and channel 
inverference can also be gauged indirectly from these statis— 
tuCs ss (6.2. Lhe umber "or [P/O watts during @ particular time 


interval). 


Pe4.1.5 File Characteristics 

Vibes tare honmalivs orecani Zed siniWays Giew are moa. 
Preceienu fOr eheir Svoregce and revriceval. Magnetic tape 
Piveswore ss eneral yw CYmat led eWwiuleba xed OF Varilapleomenctn 
CeCcordse that can be blocked.5 MInescholeeybervween, Usdtic fixed 
Ce eVvar rable we Ze records 16 (Usuallyamande on tne basitesol une 
PioLewACOLVLLY anGe une PECOr dist Zen Vara Oy ey ASmoUC peo nis 
design parameter is considered independently of the bDilockineg 
Criterion, and herce does nov enter Antosthe mathematical 
Lociitrapwon of this Optimizations problem: 

In the case of a system redesign, the existing 
Tneormar vom On blocking Lacvors (Xs); record sizes (R,), 


block sizes (R;X normal volume of files (Fe), and the 


j) 
number of alternate I/O buffers (a,,) per file per proeram 
WiLieserveras a pasis Tor the a. priori sestimauion of perform 


ance itainssand Che a posteriorisevaluavion Of sccst—erlective— 


ness of alternative system designs. 
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Lee OmOOCteMeree Utilities 

These utility programs very often impose an implicit 
constraint on the maximum practical block size, above which 
the system performance becomes significantly degraded. 
(G-2. in the case of disu sorts. the sorting efficiency 43 
Signi hrcangily reduced Lf .the block sige, exceeds one track on 
tne Gdisk, which is 3625 bytes om the §1BM/2311 disk or 7204 


bytes on the I[BM/2314 disk.) 


Weeiwe Desirable Information 
The following categories of input information are 
eensidered highly desirable but not essential for the system 


analysis. Some of them may not be easily obtainable. 


ieee eelans Lor hupurer Growth 

As far as possible, the analyst should take cogniz-— 
ance of long-range and short-range future plans for changes 
in the hardware configuration, software support, application 
programs or job scheduling. “These plans Mayeranse trom 
Bpecific snort—Lerm modilications v0.11) —det ined and ssweeping 
long-range objectives. 

Take for instance an immediate plan for the acqui- 
sition of an extra memory module. This would leave ample 
room for a major revision of the memory constraints of each 
subsystem. A ramification or goal of this decision may be 


the rescheduling of entire application subsystems, facilitated 
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by the major improvement in system performance. Other changes 
miehiuybemmorestocelizedsasuchtas the rescheduling of aroarti- 
cular job into a larger memory partition, or more sweeping, 
BUCH, 2s. ahmajor systemsereconversions 

Obviously jVnoenardsandifastirules@cankberlaia down 
on hewsmuch’ significance’ should®be* attached *torcertain antici- 
Patveds changes witherespect tolthestformudati onmot tatmacvhematical 
mode@y of *thetsystemipewhat canvandrshoutdtberstressed inethis 
regard is that any system analysis must be sufficiently 
CArCUMSPeCLIVesOny UNnUSUa le condivionssoreunanticipared changes 
that could make the’ system inoperative: These considérations 
may be the least tangible or amenable to mathematical formul- 
ation, but in more profound ways determine the actual payoff 


ana? andeed: thevoveretional) lifespan, of a .systemodesizne: 


He 42h 200 Peak Loads. 

Most Systems ere designed) wWitheenoucn "capacity to 
handle estimated periodic peak loads. The drawback of this 
approach is> the Vess=-than-fubl utilization of avaiable 
resources during non-peak periods. Other systems are designed 
Vor fanadwe the average volume mostreriicient ly. wiltir the excess 
Jeadhatrepedk pericds handled inemiicilently» ortsi phoned sori to 
be processed by another system. 

The voluneror transactions a sisuchalcace. In) pote. 
During the cycle time of a subsystem, ‘a transaction file or 


aimasver file may grow or subside: in Sizes Wheparcenn of 
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Vetrlavlon Nay navesperioddc characteristics oversaoal live, 
Weekly; montihiy, Oc a. yearlyi cycle. FPorsinstance . a tbrans— 
aC lLOnelileamayeerow sleadiily, from >, 000, records 4m) January 

to 9,000 records at year-end. At mid-year it may have an 
avcve react ZesOr (O00 records. wll slits te wSaZe i Ls 
Pebimaved, ab \Whieemid-vyear ii purcvor 000, records. (andl is 
Opximal biocking factor, iSmcalculatedsbased on this, then wats 
ij/O ,eitaciency is expected fo be preater. Than, what. at should 
be during the first halt of the year (ab. the expense, of other 
files in the subsystem), and less than optimal during the 
latter half of the year. This phenomenon is by no means 
Bisque. to fide design. But it does represent whe kinds os, 
eompromise that will have to.be made. in many instances in 
Ordere tO achieve-amoverall optimum. Other. dmportant, measures 
Oiepean loads include tie Tun Cinesper propramn. Dera ollos Vy saeul. 
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Lele Oyetem Oar kenecks 

Bottlenecks in a system may be classified into two 
Beneral categories, namely injperformance or in scheduling. 
Performance bottlenecks are considered more localized. They 
can occur in a program loop or segment that consumes the bulk 
Of the processing time or 1/0 resources Of a program or sup— 
system during its execution, such as a sort program. sSchedul- 
inpwborcleneche,eon the other band.) tend to sbe, more eencralazed: 


They can be a tight deadline which a job or a whole application 
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BSPexPeCved sto Meet, orsthe introduction of an additional 
DrOosram 1nko the production job stream of 4 tight work-shirt. 
inese considerations definitely intluence the formulation tof 
thermathematiceal modeleand thesa posteriori adaptationtofits 
SOUT ao. 

In a paper on multiprogramming systems performance 
measurement and analysis, Cantrell and Ellison [ 6] presented 
some poignant observations on the crucial effects of bottle- 
nec beyonelheppertormance xo RaMsy stem pel neyecan DCm ust ueared 
byetherfollowing examples: 

Aeiecertain trequentuly Used) program mileny anvolve 
Hour tape Liles wunnineg in a system with four tapemdevices. 
Ties wnpup sand OutpuTr Master files both, have erown Uo) oe cupy 
a litilesover ene fU1ll reel. UThe run time Jor this®@job would 
increase by a significant amount due to the extra tape rewind 
and set-up time. This problem may be avoided by introducing 


a tape wength constrain, 


iter, 

x(G + D xB) < Ls rt) 

X 

J 
where G, D and B denote the same parameters as in equation 
aia. wandelLedenotesatherieneth tof abiuairrecl (less the length 


mequiredmtor dabeliinformation) \illm aesimidar manner seone 
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Ena ute oa soya <9 JHE, (4.8) 

Ke, J 

J 
where S, C and B denote the same parameters as in equation 
See wpa no a denotes the total I/O time which must not be 
exceeded by file j. 

In trying to incorporate these two new dimensions 
of, constraints, namely length and time into the original 
memory constraints Matrix, a .mixed type constraints matrix is 


thus created: 


Sr eaten \oniias) hn 0 + 0 < My 
See Re pat f fanesh OX ata 0 + 0 <M, 
(4.9 ) 
F. 
0 Saeetowah 0 + ~1(G+DxB)+ 0 <L, 
X, 
j 
By 
0 Se Ud 0 + 0 + —(S+CxB) < T, 
Ay 


Besides introducing new variables as well as extra constraints 
into. the model, this extended matrix combines linear and non- 
idnearm.constrainte. his: mixed. type) meathiog requires more 
elaborate methods of treatment. Of course, a transformation 


of variables can be made such that a constraint’ becomes Ianear. 


1 


i 
J _____ > —— (4,10) 
By (G+DxB) L. 


i) 1 A's 
7 _ , . 
= - | 
7 = 
- 
(Dat) 
. 
ve cr awes ion “ie a6 aah =) Sue i ae 
® ‘ 7 
> 7 
‘oe 
av ‘ 5.3 aniven 4 7 gos. add | eeHOReD * bin > 
-f ote? ew babe 
' a 
mi fl a i e =} 3 4 +- 1 : 7 
ey 
\ if I anes Vigne s Pees a | 
: . -_ 
} a} ' » 7 A ) - gannoe. wee 
7 
— a 
: - once nant! 
- > na és "' 4 ‘Oy + 
0 ies , rt v1 sam 7 
} af A A 
e - iz 4 
SS a ij- & 
. a ; y 7a "\4 ae + - *. —_ 
A 
uo q - 7 
wWiHiwet L197 « : jas US. 20 Palen 
nnd? . (ens aa ai Mm waite ae | ings is alae: > e aie 
; StS" OCTSOGTS 8957 Et 


yy 


wl seme ab : lc aha WwW 


o 
tintoriiaas » sk 


a 


— 42 


But this does not alleviate the problem of solving for a 
unique xX, for as particulars filemoccuringey inmidirfrerent. types 
of constraints. Furthermore, a composite objective function 
would need to be defined in accord with the composite cons- 
traints Matrix.= The) different: components of®any such overall 
measurement of system performance would have to be weighted 
enipimiiical lye wine this veeneral®@ area LA. OSilverpetsal [29 | 
have proposed a weighting algorithm to aid the systematic 
allocation of system resources. 

In our present context, the above-mentioned refine- 
ments would create an excessively elaborate and unwieldy 
mathematical model, both from the points of view of system 
analysis and Tanding- availaples algorathms Tor itse* solutions 
The payoff for this elaboration of the problem is not immed- 
jiately apparent. Indeed, the complexity may deter the manage- 
ment or analysts from undertaking this project. 

However, one can take a crudely empirical approach 
to achieve the desired results. In substance, one can arti- 
ficially bias’ the’ parameters! of the affected Pileshin’ the 
model. The best predicted solution is then adopted. It should 
be re-emphasized here that the overriding objective is to 
achieve useful results with the minimum necessary effort and 
a sound methodology. With due regard to its effects on the 
physical system, a mathematical model can be distorted in 


order to simplify its algorithm for solution. 
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4.4.2.4 Pile Activities 

Application files normally vary in size and pattern 
of activity, often in regular cycles, The relative amount 
and manner by which the file sizes vary affect the coefficients 
el hvnesobyecrivestunctionedinectly .esHenceytany meaningful 
estimationvortashighlyqvariabbersfiidle!size overtasperiod of 
time would need to take into consideration the pattern of 
Variation. 

Doesythe Hitersizesl huctuatechwithas Uniform or 
normal distribution, or skewed towards a short and rapid 
growth followed by a long and slow decline, or vice versa ? 
Does the mean size coincide with the median size, and which 
One wermore representativenveaWwhatsisnthesvariancesotca tite 
Seer over saeperiod Ole time, essay, auyecotar eWheleicethesnatvern 
OUSaceLVvloy inethne Masver Wi osm es Orentne CransacuLons uni 
pornivediscpibuted, GorseiUsteredsaroundsaremadsapercentage 
oF Nich actavatyeaccountseca wnt hers tat vermicasesmither majority 
of the master records can be read/written 'on the fly', that 
is,iwithoutenecessarily waiting, for theytapermotions tor come 
tO ancompletenstopiin thei IBG.. By itself, thise consideration 
doesanotedirectiyl affect the mathematical solutions fer the 
epiimelwolockingy factors, SHowevereether tiles vee coeiiiciene 
otbavlow ectivity fildeacouddsbe artificial lyprediced wortake 
advantage of this fact, (under the assumption thatea relative 
inerease in 1/0 time resulting from a decrease in the blocking 


facthorafor a low activity files would» bes morey Thanscompensated 
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bye ches increase lom overall o efficiency of the subsystem.) 


Wao. 5) Stabistical 16/0 Brror Rates 
 wtatisticalvinrormationvon I/O error rates is norm— 
ally collected perarile volumeranditape "device 4d This Tinform-— 
ation can be used’ to gauge the relative significance of ‘the 
number VOrYer ror Iblochimetricstfor@alnantbicularetape Piles 
Theoretically, the larger the block, the higher is 
Drespyobabi bity of poccurrence fof Panli/@ (error dinkthat tb lock: 
However, the total number of blocks in the file is reduced. 
im any e@case,;tatreasonablyewell mainteinedatape installarien 
is not likely to encounter a situation where the number of 
b/O@error blocks ipertiii he sbecomesEsieniit icant senowgh*torafiees 


the formulatvionvorathe objeclivelfunctionrwith error (terse 


Mas sOeBselechioneGriterta tor Blocking Factors 

In order to make a meaningful comparison of the 
performance improvements that can be expected from a mathe- 
matical solution, the analyst would need to know the criteria 
byuwhiiennthe tblecking factors werenehosecnmnethespast. taWas 
tPiemapproach purely ad hoc’ * As is cften dikely, each suo- 
system was probably given the blocking factors that represent 
the best guess of a systems analyst. In a situation with 
system constraints, this is seldom the optimal, arrangement, 
and a sensitivity analysis of the actual and the optimal 
arrangements would indicate the amount of deviation in system 


performance. 
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CHAPTER V 


EXTENSION OF PROBLEM 


be 8Direct Access Sucrage Devices 

The foregoing chapters have identified and treated 
a problem in optimization of computer system resources in a 
relatively straightforward and predictable frame of reference. 
The senvaironment considered was atnon-overlapped)))singtetbaven 
Precessineswitneuaperriles .“hingresality, @rew Catatprocessing 
SvevoMseliuUr such atbsimplilitedsdéscription.. With the addition 
of direct access files and the introduction of more powerful 
Operating systems, the question arises as to whether the 
effects of varying-the blocking factors on’disk and tape’ files 
are mutually unrelated and additive, or are related in some 
non-linear fashion. The remainder of this chapter attempts 


LOmcOagresseunisyproplem. 


Do) Disk tile Organization 
Disk filestaret generally organized in one of the 
four methods: sequential, partitioned, indexed sequential; 


or direct. 


Bee. Sequentialborcanizapion 
Invatsequential filey records aresvorcenized solely 
Siathe Dae lector thelr successive physical locations in whe 
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sequence according to their control keys as well as in physical 
sequence. Uhe records are usually read or updated in the same 


order in which they appear. 


heecwes Marta taonede Organization 

A partitioned Tile is one that is divided into séev-— 
eral members. Each member has a unique name. Members may be 
called by name for processing. The records within the members 
are organized sequentially and are retrieved or stored succes- 


sively according tonphysical sequence. 


bie oe) Indexed) Sequential, Oreani zation 

An indexed sequential file is similar to a sequen- 
toa lpi itemineg thaterapids sequential, processingy iss possible: 
indexed sequential organization, however, by reference.to 
indicessacsociated=with, thesetide ~makessitealso possible To 


Guvekily locate individualerecords 1c random processing: 


5.2.4 Direct) Organization 

A Tales oreanized ineamdivectemannéreis characterized 
by some predictable relationship between the key of a record 
ANG the. addresses OmviopemecOrd One the Girecr access storage 
device (DASDD)Rr Mhts welatdonship is normally estapli shed by 
a randomiza nese comichm that converts a record key sintowa 


Unigues DASD address. 
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Doe Disk Scorage Orcani zation 

A disk device essentially consists of a number of 
mecOragbne tracks per disk suriaces, Hach (set of tracks thas 
can be accessed by the read/write heads by the same positioning 


of the movable access arms forms a cylinder. 


Dota rack hoOrmats 

Information is recorded oni all devices an <a forma. 
which is prescribed by the control unit and which is identical 
forrall devices. Each track contains certain 'non-data'’ infor— 
mation as well as data information. The non-data information 
generally consists of track alignment index points, physical 
Pape, eblockeaddresses=and=recordmaddress. Ihe Cara-records 
may be formatted with external keys (count-key-data) or with- 
Sutsexterna keys (count=dara) oo hieure “it liustravess the 


track=Lormats for one: partacuban-manuracvurer. 


Doe he Cord srormlatis 
Logical recordsimay, bescreanized win one of give 


formats (ac illustrated in fieures> er 


—- Fixed unblocked 

- Fixed blocked 

- Variable unblocked 
- Variable blocked 

- Undefined 
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Index 
Point 


Track Descriptor 
Record (RO) Data Record (R1) Data Record (Rn) 


pie LS Cac RN pie sd ced pei Oj Vea shal oy 
Count! | Data Count Data Count Data 
HA GIAIG 
eg [rae [Rertlc [pee | fel Cla SMS SM cla lelaialic Area 


A. Count-data format 


Index Point 


o) 


Home Address 
Track Descriptor Record (RO) 
gy Deta Record (R21) i Data Record (Rn) 
Feeg Data Count Key Data S Count Key Deta 
vo feaje [ese Area =.= Area Ic Area lc be Area Ls Area . Ss | s[alc ees | Area Ic Area 
Grants ae ravcrard 


ere Check | in t Cyclic Check 
Las Data Length 


Head Number 


Cylinder Number Key Length 


Record Number 


Flag 
[sa coe ei Head Number Identifier 


Cylinder Number 


Flag 
Address Marker 


B,. Count-key-data format 


Figure 5.1 DASD track formats. (IBM) 
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Variable, Unblocked 
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Count Key Data 
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Dries Eifects, of, Blocking onyDASD tiilegsS ystems: 


Asesone) “can »observe tintin cures slyephysicaleeans 
analogous to the inter-block gaps on magnetic tape files also 
exist between records on DASD files. It would thus seem 
natural to apply the same methods of treatment as for tape 
files. In fact, however, the operational characteristics of 
DASDe are! sipnificantly different: fromethase of ,taperdevices. 


The two basic device characteristics which distinguish DASD 


- eee : 


aqeq Dpotayita) £12, oust?! at veseade agp! ane oh 
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from tape devices in blocking considerations are their constant 


speed of rotation and their fixed track length. 


B41 Eifeectsion Pile: Compaction 

THe primary reason-eLor plOCkKane recCordseitimLASDaismaL es 
is to pack the file storage more efficiently and hence reduce 
Slorage space vequirements. | With blocked recordss | raps exise 
between blocks of records rather than between each individual 
bogical record.) However, it should) be noted from fisures 5.4 
that the file space requirement does not decrease monotonically 
With ancreasing plockingsractors.— lhe us) dle voO= Thee lace 
that a block cannot be spread over two Prato ee Maken, ullese 
Cemulain blockinewiactorse ss the last, block ton ar track overt Lows 
one track and the entire block is pushed onto another track. 
This resulte: in some wastage per traci wand ie ebecomes smore 
significant as blocks become larger. These relative increases 
in the amount of storage space wastage are represented by the 
local maxima in figure: 5.3. .> The same typestor characteristic 
curve applies to any track-formatted DASD (i.e. disk or drum) 


which does’ not allow a block to be split over two tracks. 


12 Some DASD are equipped with a feature for splitting a 
block: over two tracks. In such a scase, the characteristic 
Gurve will resemble those of tape devices, i.e., there are 

no local maxima due to track boundaries. However, they still 
exist due to cylinder boundaries, though fewer than before, 
INvuicively. 


-_ 


- ome - a ; 


nia 


78 ite = ar 40 ain _ We, we 
; he eo ens 

: do eg ganeey ae 2 5°: 

Aree See ee at oe eee 

i pee wud . Sir ae 

: : ' a ay * 


Pe 
a “joisg ohm Me ope heated 


salt? I val wr Ses ty “37 19° .e Rot. wena “<3 apt 


- 
* 


: ‘ a =e 
=e A 7 , ort 0 he ba +E navi © bat yroto ale éat ‘Neen 0! @t 
4 ’ a a a a ha 
Pe y 
4e.) 
‘ ? ~ <= Ga 14 ¢ are | L 133 
fa4é Dt oo. 
a 4 7 
o> ‘ f 5 an j ~ * 7 + ® sive ri 
* * a 
yy ci ) } Tek. | : 
¢ ic a wT 
‘ 4 : 
. coll i> ig ears 
ay i . 
, ‘ i 12) - nar 7 ~~ 
x ; a 
c * | nae =I ‘9 an , : aa ia reg 4 ‘ 7 
* * 
+ ‘ 
“is i 7 oy = ee ‘ a a 4 as arns srs ars rece 
* : Sd _ . : 
all : - 
cH mibad $4 Cine 2 whtTus SapAey se Hf 8? easy 
‘ 4 4 . i® ’ 
2 >» © 7 ra = =p ay wi 
WRAMLN! GYI,S leet. |. ese! Sincet - cuenta. ce wyeet 
c oa . . 4 * -* 7 ; : 
‘ eo ¥ 
> : 7 : 
+ Z EA at - et 4 «) e@hi;% rs) re n Ipare 
- ; J / 
oy P ; - - : 
LIfe ’ 70 +4 ‘ 5 | i wy Ay ate | ee i 
& - ae , ‘ > . 7 on 
im th 75° ; ri Lea 2 ! = ee +) a oe oS q 5 ® e an wake ae “a 
7 ; sd cf bs A i‘ 
.” : a» . a r : 7 
= = - 
, abe a W TU é ’ dl i) — som eK. Weld 
7 a » "4 > . : zal : 
i Kr 4 : os : nae 


oe a 


aise) 


BLOCKING FACTORS 


| 
| 


| 


| 


| tH 


SNIVUL 40 “ON 


FIG. 5.3 NO.OF TRACKS REQUIRED PER FILE AS FUNCTION OF BF. 


(IBM/2314 DISK DEVICE, WITH KEYS) 
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5.4. ete cecte on £/04.Time 

AnOtherereasonafor biockingsds thathitimayasave 
titecuyeliurecordss arerprocessedyconsecutive ly, yiherentstonly 
one rotational delay before reading or writing a block of 
records... yHoweven, if | records: aresnotsprocessedyconsecutively, 
blocking may be a disadvantage, since it takes longer to 
transfer the entire block rather than the single record to be 
processed, 

Asmore’ compelling) reason for-biocking records on 
disk devices with movable heads is that the file compaction 
reduces the total number of cylinders required by the file. 
iiecheslDASDetides 1 sSwassignéed agerouple?r convaicuous. cyiinders., 
tha seereductitoneresuitve tingapsmaidersaveracesnumbereot cylinders 
traversed by the read/write heads. 

im ancdirectily,organized disk) file) Awhere (Ghesgrnandom 
input transactions require a short response time, the trade- 
eff between an increase in block transfer time and a reduction 
in the average seek time depends on several related factors, 
namelyeathesdistriputton of} ¢highsactivityarecoras tnethestite 
ahdethespivei cal —sscatterconethevitlesimntognon—conctiguous 
cylinders. Here one must distinguish between disk devices 
with movable heads and drum devices (or fixed-head-per-track 
disk devices). In the former; the block transfer timesis less 
Significant than thevaveragercseek ime, sof auratic cial itoms . 


In the latter, however, the seek time is negligible (at elec- 


tronic speeds). 
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DeD Leolbationror orropten 

From the foregoing discussion, it should become 
obvious that not all forms of processing or methods of DASD 
file organization lend themselves to the straightforward kind 
of mathematical modelling that was applied to tape file systems. 
Por instances an-atdirectily*organized@drim) (or Pixedehead disk) 
file that is processed randomly, the optimization problem 
hinges on a space-time trade-off involving a qualitative 
judgement which is beyond the objective of this optimization 
Problem, birenecoOmmMinimizest/ONGimetor steraceespace™ plbisaisc 
nov’to say that the method of treatment cannot be further 
evaboravedttoedealewithethis particularssituationyeabut for 
thespuUucgposes Off thisaidissertationsrthisecacescan besisolated 
and removed from consideration. 

We shall therefore only consider DASD files on disk 
devices with movable heads. The working assumption made here 
is that the reduction of the average seek time more than com- 
pensates for the increase in block transfer time when records 
are blocked,optimally-*""This approach would thus apply also 
tomeandem processing Cf UWirecilysorranimecnl ied. 


Poreréasons of functional modularity and taneiytical 
Simplicity. ~the disk Tile subsystem will pe considered Wepa— 
rately from tape file subsystems. Their separate mathematical 
models may lead to sub-optimal solutions. But the functional 
disparities between the two types of subsystems often necess- 
itate different requirements for modification and reorganiza- 
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5.6 Methodology of Treatment 


5.6.1 Formulation of ObjectivesFrunction 
The*objective function for disk® file subsyspemesis 
formulated differently irom that of tape =filessubsystems. 
tne primary .~Obj/ective NHeresis sto minimizerthe? total number 
Gh tracks requiredsby avdiskitilet= Because,or" thet fixederrack 
length, the minimization of the total ‘number ofr disk record 
blocks does notenecessarily imply aominimization’ of the rotal 
number of tracks@requiredre this fact is illustrated by the 
existence of locad maxima on the graph of [igure 5.3: 
Theatvotalinumbervor tracks tosbesiminimized for n 
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#. The fractional parts are truncated aiter) the operation. 
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where Py pie X, have the same meaning as in equation 4.6, 


J 
meGsulotal noe; ofttracks requinedtfor asdiskiiides 
Y, =eaNomsofeblocks Sper trackr for (dieksfi tess 
UO @> Track capacity sinenomofecharacters 
MS =ahonrothcharacieps sin Glas bab heckfon hrcack 
W, = No. of characters per block (except last block) 
= Key Tenet an nO wOmecnareacters 1Ore tates, 
B, =sNo niet acharacters penpbilock iiior fides) 
C = Device constant (C=0 if eee O20, ise K #0) 
P. =sdevyrecenconstantecallows for positional deviations) 
Q = Device constant for non-data information in 


no. of characters 


5.6.2 Formulation of Memory Constraints Matrix 

The memory constraints matrix for a DASD file sub- 
system has the same form as in formula i 5 Bub te ot nouenivd & 
has the same form, more care should be taken in setting it up 
in conjunction with the performance criteria which are not 
quantified into the mathematical model, such as channel con- 


figuration, DASD file organization and access methods. 
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5.6.3 Diserete Optimization Methods 

Because of the nature of the function as illustrated 
Dy Une loca lemaximea in Bie. =5 23,0 che sopcimizatvon problem is 
best solved by using integer programming methods. [90, 91, 93] 
Thiisnelass O1 aleorithms is peyond the scope Of this thesis, 
and hence will not be discussed beyond listing some key refer- 
ences. 

A possible alternative approach is to approximate 
tiesaiscrete funetion with a polynomial and then proceed. to 


Opuam~Zerrne Continuous function. 
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CHAPTER VI 


IMPLEMENTATION 


6.1 Research Recommendations and Managerial Action 

The subject of operations research implementation 
has been treated in depth by Huysmans [52]. In his book, 
Huysmans stresses the importance of the involvement of manage- 
mMentway alin levels 2n order to Traciiitate a smooth sand success— 
ful implementation of the OR recommendations. The failure to 
secure top management agreement on basic objectives and scope 
Of the project at the initial planning staze may jeopardize 
Lhe eccepvance Of the recommendations at the amplementartion 
phase. The possible behavioral reaction of the manager to 
the research recommendations as a function of managerial 


imoderstvanding, can be clasesifiled vinwone or Tour caterories: 


- rational rejection 
- resistance 
- acceptance 


- implementation 


6.1.1 Rational Rejection 

The manager understands the research recommendations, 
but rejects them on rational grounds. The research may be 
imeadedquate, irrelevant, or nov superior (oO presenl practice 


GLranertaonm coct considered). Irrelevance of the research 
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Giver takecmeneesorm of an improper consideration of a sub— 
system of Une potabasysteme.s Thatvis spthe wreongeconstraints 
bave wbecwesinel ed OUT Al oreoptuami zation, 

Rational mejectiloneis related. to, the problem of 
pre) ectheselectiLon, that iss, the choice of a subsystem to’ be 
considered in the research. It should be realized that the 
evaluation of potential implementation problems should be an 


important Gererminant Of sprotect mse lection, 


6.1.2 Resistance 

the Manager srejecus the research recommendataome, 
and his way of thinking about the research problem deviates 
Considerablyahronethae Chethesonalverls, Rejection in this 
case may very well be caused by a lack of managerial under- 
sLandine®. Sali evnewanaly sis miseadequate, relevant, and superior 
to present management practice, there exists an implementation 
problem. The manager may fail to recognize the advantages of 
the research recommendation for various reasons. He may not 
Have paidusuffictentdattenbionstosthesrescarch propesad;eike 
may incorrectly evaluate the proposal's return, cost, or tran- 
sition-cost; his»personal, goals may conflict with and take 
precedence over organizational objectives; or he may have 
insufficiently resolved the inconsistencies in his perception 
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The manager adopts the research proposal, but does 
not understand it. An implementation problem exists to the 
extent that instability is introduced by mere acceptance of 


the research recommendation. 


6.1.4 Implementation 
The manager understands the research proposal and 
SOOD CS. Li. meLoYels USseiulL LOrUlymece Unis —Cabverory anv —GWwo 


sub-categories: 


- sustained implementation 


—- autonomous implementation 


Gay PME @ Sustained ImpVvementation 

Managerial understanding of the research recommen- 
dations» is integral, that-is, the manager™has good overall 
undersvandine of the critical factors, underlying The r,econmmen— 
dabions, but this understanding as acquired au the cost of 
continued involvement of the analyst. The manager realizes 
that changes in environmental parameters may require changes 
in his actions. He is able to identify the change points, 
but proper reaction to the changes is only assured if the 


analyst stands ready to provide his continuous support. 


4 Ou 2 ’ 4 ; Hi a ">t i Re 7 = mye: sa 
: at. 
Ser 37.6% 7 grates aris lfetesnelae:- he, ot & 
= a : ; } § 
v 

_ - - - . 
als ite Lael vaersaae ats 


yy 


s jos 7 
_ “vs 
- 2 FA) 
> 
: ‘ . 
7 
3 
pic, 54 - 
= > a | 
% 
— ad 
» - 
+7 sr =" Se - 
- — ee in 
7 - 
" * _— - ® - 7 
“ j = ped © a s » f a ten nm a ? 
-.. keene 
aie eA rae > -¢ ad inh eee 
¢ —? oy 
ae | é en = « ’ r 
‘ "? 5 
- 
me: Bat ae mer j o& Eir ss *netim Vv Loven hour? Se 
' ’ ts 
>» = 
: a e° 
wigeeean ee tae Vea a4 s Vee aah | yer ares ¥ > SeRTAKy. 


os : = 
; as 


: z- s 7 
saihety aGante od? inet ot etawex lel sanifume 
= - a _ 7 - 


- a ) 7 - 
@ca~Ag taviiags (dem 2) eopnen’4 cat 62 neervee comiea 


— ~ 


. a _ : ; _ soflt) - , 


7 
Fags, aieabties ott. 
_ ‘ : 
. i 


- 60 


6.1.4.2 Autonomous Implementation 

Managerial understanding of the research is explicit 
and complete. No continued support of the analyst is required. 
AdopuLon occurs because athe analyst's proposal is consistent 
With the managers penspective, and solves .a problem forshim. 
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6.2 System Performance Monitoring 
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its performance is monitored .and evaluated recularliy.. Atwany: 
suaccuof thems ys veme bi femcycle ~wehesnecdsback from tone phase 
may necessitate Fae mee oteay Oe preceding phases... The .moni- 
tored performance of a subsystem thus ultimately determines 


Giewvaltdicy and reliability of this sORsapproacn. 
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CHAPTER VII 


CONCLUDING REMARKS 


bet Economics Payor slvandeintanr ible (Benefits 

The foregoing chapters have described an operations 
research approach tTomsolve a) problem tof optimal resource 
allocation in tape and disk file systems. Analysts who are 
convinced of the immediate payoffs in increase of throughput 
and reduction of file storage requirement still face the task 
of convincing the management of the advantages of this mathe- 
mavecal Optimization approach, ~=lhils “approach “te problem ssolv— 
ing will have to be sold to the management in terms of manage- 
Nene's Pramesof reference’ —— namely; "systempobjectives;s pro- 
ject development costs and economic payoffs. 

im this bight f’theerile compaction “advantage may =pe 


transiatedrinto the follliowing managenient “considerations: 


- reduce tape and disk acquisition costs 


- improve job scheduling 


Simi Vandy e “the *reductiom-0bSi/OPlime andFovernéad may wbeetrans— 


Wecved@intomtne following Management berms: 


— increase throughput (or reduce computer time cost) 


= reduce bottlenecks “in Job “scheduling 
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What may be equally valid and perhaps have more pervasive 


impact are the intangible benefits: 


- greater systems awareness 
- improved resource allocation 


- better management control 


(1-2 Applicability lin Fucure Technology 


The costs of memories .at’all levels are continually 
decreasing as technology advances. Already, it is reaching 
a point where the price/performance ratio at each level is 
becoming competitive with that of the next lower level. Some 
installations are now using large scale main memory for data 
which used to reside on DASD several years ago. Disk and 
drum devices have also evolved to the stage of sophistication 
where the price/performance ratio makes them competitive with 
the old magnetic tape technology, which itself is reaching 
new heights of speed and economy. 

However, with £Uull recogenitionvor this rapid vech— 
nological obsolescence that characterizes the et gas aby qh Riis ohio 
the upgrading and diversification of hierarchical memories 
Will not invalidate this mathematical optimization approacn 
to computer systems resource allocation. In fact, the present 
trend of increasing complexity in computer systems can only 


lead to a greater reliance on such analytical techniques. 
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LocaL OPTIMIZATION BY LAGRANGE MULTIPLIER METHOD 
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APPENDIX A 


LOCAL OPTIMIZATION BY LAGRANGE MULTIPLIER METHOD 


Let us consider the following example of a Master Update 


program: 


F,=50,000 L.Recs. 
Ry- 150 bytes /irec. 
ay a2 alt. buriters 


Trans 
action 
tape 


Input 
Master 
tape 


Pz-107000 LeReCS. 
Rz=80 BYTGsy rece 
diz = 2 AleewoUa nets 


Master NO Gaia Dit heise 

Update available 

program M = 10,000 bytes 

i 
Fo=50,000 L.Recs. Report F,=20,000 LW RECS 
Reais O0a bytes] rec. Giske Ry, =70 bytes/rec. 
ag=2 al tebita ers aye aiteiwlidht ers 


bisa Ang application, progham,ustnewiounstabe ules. 


THe Opjective TUNCLLON a GWwita Ney Gwe GoGtLOheGrOllta mrs 


ae ee (i) 


The memory constraint equation is 


2 x (I50KXy re Lo On AG 80xK, = 70xX,)) See 0 U0 mee tte) 
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Taking partial derivatives of (i) and (ii) with respect to Xie 


(j=1,2,3,4), we obtain the simultaneous equations 


- 2 + 300 X= 0 pe 

x2 Cilio 
ih 

. os + 300 X= 0 
~ 

= fe FIG 8G & 
Xz 

- oy 2 a 
Xi 

Solving the equations (ii) and Ciii) simultaneously, we obtain 


the mathematically optimal combination of blocking factors 


X,=12.1, X,=12.1, ey) VG CHa. 


which we can round off heuristically to give the physically 


optimal combination (occupying exactly 10,000 bytes) 


Xy=Xo=12, Xz=7, X,=12 


From this we can determine the total 7/0 time for the four 


files 


Te=el58,,58 +. 158-58 ete he/ Sete she le SePCOnGs 
% 
= 378.86 seconds 
- Example using the formula for !BM/2401 mod 2 (800 bpi) 


Q9-track tape device. 
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Sensitivity Analysis 


If the BF's were arbitrarily assigned, say, 


XpRX5=Xz=Xy= iki. 


(with a slack of 100 bytes to spare), the total I/O time 


would be 


TOL=G] GIF61 ET GIs. 6 PEA), G50 SHRO2 Seconds 
= So UssSeseconds 
which differs from the optimal T by 
Ae Eh 0T Wy £7) / T= 0.8% 
If ton ithe ‘other hand the total core memory avail- 
able for tape buffers were equally divided among the four 
he Pees 


fites@Mthis) assigning 9Xj)=xX, =8, 9X.) -15, Xs aie 


(occupying exactly 10,000 bytes), we would obtain 


4 175.25 + 175.25 + 18.69 + 31.38 seconds 


400.57 seconds 
Whichmdiesersesigniticantly trom the optimal le by 


Ag 4=)66T vs), vl? 4) Ore dk 


The ditferent degrees of file compaction for the 
ahbovescasesucan also be calculated using a similar formula. 
It should be noted in passing that the percentage differences 
remain the same for the same relative difference in file 


sizes. 
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std 10% aobsoaq@os ski lo zeorgsh aaa | 
-nsniead aah ae ie eich st ca eae 


C PROGRAM TO SOLVE FOR OPTIMAL TAPE BLOCKING FACTORS 
C BY LAGRANGE MULTIPLIER METHOD. 
C 
C FORMAT FOR INPUT DATA: 
C CARD 1 
C COES? 1—3 NO. OF OBJECTIVE FUNCTIONS (INTEGER) 
C CARD 2 ONE CONSTRAINT PER OBJECTIVE FUNCTION 
Cc COMST 1-2 3 NG. tO MIAPE FILES PUINTEGER) 
Cc COLS 4-13 MEMORY AVAILABLE FOR TAPE BUFFERS 
C CARD 3 RECGRD SEZES 
C 10 COLS PER FIELD, 8 FIELDS PER CARD {DECIMAL} 
C CARD 4 Heal ee RASS 
C POR COUCST PERG TEED, cunt leLDSe Pen GARD SUDECI MAL)» 
C SUCCEEDING CARDS REPEAT FORMAT FOR CARDS 2 TO 4. 
C 
DIMENSION X(9)5 RUG), FL9), TI9), CON(9) 
READ(5,101) M 
WRITE (G6, 100) 
DO 200 J=1l2H 
READ(5,101) Ny C 
READ (5,102) (R(I),I=1,N) - 
READ(59102) (FC(I),I=12N) 
C BUEEEReMCMURY/ Zan COUR S2 ALTERNATE HUEP ERS SPOR] PAP EEET LE 
CK=C/2. 
S=0. 


DO 130 I=1;N 
T(LTI=SQRT(RCI)D*FC1)) 
130 S=T(1)+S 
W=5/CX 
DO 140 I=1,N 
140 CON(T)I=HR(1)*2.6 
(RITE(6,103) {(CON(L),I=LloN),C 
WRITE (6,104) 
DO 150 J=l,N 
X(T )=SORTOFCI) I / CSORT(ROL) } €W) - 
150 WRITE(6¢105) LsX(l)gTsRi td eTsFCI) - 
b= Ve 
DO 1380 K=l9N 
180 Z=Z+F(K) /X(K) 
WRITE(&,106) Z 
200 CONTINUE 
190 FORMAT('1*) 
1iG1 FORMAT (13,F10.0) 
102 FORMAT(3F10.Q) 
103 FORMAT(//,*GMEMORY CONSTRAINT EQUATION 3145//316F8.0) 
104 FORMAT IS OCPTINAL Bek «UU Se ys Peis RECURDUSIZES<..3 
SAT eRe Sop aes ee | 
105 FORMAT OC! XU si2e') ="sFrGertet2ls' Rie l2s')) =" 765.0, 
M4, PEL sleet d he ere n 
106 FORMAT(C'OMINIMUM VALUE OF OBJECTIVE FUNCTION 2 ="'sF8.0! 
STOP 
END 
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MEMCRY CONSTRAINT EQUATION 


300. 300. 160. 140. 10000. 
UPTINAL SE. Fo 7S: RECORD SIZES: 
X{ 1} == P2208 R { 1) “= TUS 
X{ 2) = 12.08 R{ 2) = 150. 
X( 3) = 7.40 R{ 3) = 8Q. 
X€ 4) = 11.19 R{ 4) = TO. 
MINIMUM VALUE OF OBJECTIVE FUNCTION 
MEMORY CONSTRAINT EQUATION 3° 

300. 200. 12000. 
OPTIMAL 8.F.*S5 RECORD SIZES: 
X{ 1} = 31.79 R( 1) = 150. 
X€ 2) = 12.31 R( 2) = 100. 
MINIMUM VALUE OF OBJECTIVE FUNCTION 
MEMCRY CONSTRAINT EQUATION 5 

200. 200. 8000. 
OPTIMAL B&F. *S3 RECORDASIZES - 
Ate) ese oF R{ 1) -= 100. 
X€ 2) = 9.61 R{ 2) = 100. 
MINIMUM VALUE OF OBJECTIVE FUNCTION 
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F( 1) = 50000. 
F{ 2) -= 50000. 
F{ 3) = 10000. 
Fl 4) = 20000. 
= 11415. 
EVLETSTZ ES. 

F{ 1) = 50000. 
hig j}2= 5000. 
= 1979.6 

PILE SPES- 

fil Lje= 5000. 
F({ 2) = 500. 


= 217. 
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APPENDIX B 


GLOBAL OPTIMIZATION BY GRADIENT PROJECTION METHOD 


Let us consider the following subsystem of 3 programs 


and 6 tape files: 


Fz=10,000 F,=50,000 
R3=80 R ,=150 
Oaee aq 
(ae | 
j 
be eS Gettin, yng Mery pt A 67200 


ff Re=100 


Program A 


lFreq=l/day 


F,=20,000 


Bre ipen A subsystem of 3 interrelated programs and 6 tape 
G1 Less, 


a * 
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The objective function ts 


z = 24.25 (50000) + 24.25 (50000) + 22 (10000) 
Xy Xo X3 


+ 22 (20000) + 5.5 (5000) + 1 (500) 
Xy 


Xs X¢ 
The memory constraints matrix is 
X} XQ X3 Xy XS X6 Mj 
300 300 160 140 0 OGes LO00D 
300 0 0: 0 200 <2 ODD 
0 0 0 0 200 200m =s 8000 
1 0 0 0 0 WP eS 0 
0 ul 0 0 0 nee 0 
0 0 1 0 0 i) eS 0 
0 0 0 1 0 0 2 0 
0 0 0 0 a Oe 0 
0 0 0 0 0: i 0 


With the use of a nonlinear program, the optimal solution 


for the tape blocking factors was found to be 
Kaho eelicenl NS = Tol, att = eLOs Cw AG = Oo e, EAP om. 
which could be rounded off heuristically to the physical 


optimum 


XyeXo= 12, X32 7, X42 12, Xe= 36, Xgz 4 


Bid wgk vS.2o sek .Oe0L wah yO.0 meh (SS0 agit= x 


Festa <6 a-eS eat SAID: Tie Mena bine A 


ne rel Go nat 
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The buffer memory used, tape lengths, and I/0 times per run 


of each file were computed to be 


is ia 
394.696 394.696 76.2827 104.55 21.489 <— 
BIOL 7.00) DAO. rode oli Gel opel oo acae eee! 


IOTIME (SEC.) 
LENGTH (FEET) 
MEMORY USED 


36005536000 11200 1680. 7200 £00 
T2000 2 


The minimum value of the objective function, 1.6.5 stnemtotal 
number of tape blocks processed in one run cycle of the sub- 
system, (one month in this example,) ts 2z = 271067. 
Over one run cycle, the I/O time per file becomes 
TOTAL penal 
LOT 1ME GCSEC .)e=) 9513560574690 al Gira ec 23,0,0 411d to oO 
It is immediately obvious that any significant percentage 
reduction in the I/O times for files 1 and 2 would improve 
the system performance much more substantially than if the 
Same percentage reduction were effected on files 3 to 6. 
In general, large tape files that are run frequently 


should be given large blocking factors. 


Sensitivity Analysis 3: 


If the tape blocking factors for program A were assigned 
arbitrarily, say, X4=X9=Xz=X,= 11, 
and those for programs B and C were subsequently assigned, say, 


Xp =X¢ = 207 then we obtain the results 


1b, SPOOR AG Ta 
IOTIME (SEC. ) h07.853 807.850 58.2395 109.815" (2 oat oeee 


LENGTH (FEET) = 637.311 637.311 91.0038 171.591 39.6094<— 
7D 760 1540 4000 4000 
MEMORY USED 540 0mm O0mme. Bryan 


2.53483 
7 
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Over one run cycle of the subsystem, the I/O time per file 


becomes 


TOTAL an SEE 
LOTTUES(USECS))=9890.U2 949990742 ~ 1283.27 “24157 8oee 39.445 


Percent deviations from the optimum can then be computed from 


the formula 


q 
Daf Oi tank Gone, 
A' = rd Wot 74 
e fj Tj (opt) 


where f; and ur denote the same as defined earlter. 


Thus, 


A= AM; aa 
PERCENT TIME DIFF FROM OPTIMAL = 3.33 3.33 723.7 5.03 18 


PERCENT “LENG DIFF FROM "OPTIMAL = 79°933 ©9333. 101207 5.03 18 
PERCHEALILOTAL TIME DIEE FROM OFTINALS = 164 _57 al 
PERCENT TOTAL “LENG DIP *PROM OPT IMAI, = Fi 39 1 


If the buffer memory for program A were divided equally 
among its four tape files, and then files 5 and 6 for programs 


B and C were readjusted arbitrarily, assigning 
XpsXo= SPenXs= 5, PA sn2 07 r= 15:9 Aer pee ty 


then we obtain the results 


TOT08S 
N67 .05 7. 8675 05 J et Oe. 9 tee a oro eerie 
729 COL Oe Leo oO Oe mt et eee os) Sehr) Menara os 
2400 2400 2400 2800 7800 200 7 


LOTIHE (SEC.) 
LENGTH (FEET) 
MEMORY USED 


[ye 


ZOnwoD 
Z = 340997 
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Over one run cycle of the subsystem, the I/O time per file 
becomes 19.035 


DOD TES Cle Coeyeeaenl a 26h eee 326.1 1096708 91790,60 feos 
TOTAL 


Percent deviations from the optimum are computed to be 


we 


hee i Leyte oa 
PERCENT TIME DIFF FROM.OPTIMAL. = 18.3:.18.3 . 34.7 .022-1.</ 


PERCENT LENG DIFF FROM OPTIMAL jie 1a 2 lee OD acat 
PERCENT TOTAL TIME DIFF FROM OPTIMAL = 10.4 pi nla 
PERCENT TOTAL LENG DIFF FROM OPTIMAL = 10.8 wo bah Ss 8 Ie) 


A brief look at the sample program output in appendices 
A and B will show a distinct difference between using local 
and global optimization. Local optimization for program B 
yields i Ciba Th Skee Ms Sih, 
whereds: global@opimmizationeyleldst #xy=e1202u parX cS oe 2k 
Thus the two approaches are not always compatible. 

It should be emphasized that these percentage reductions 
pertain to I/O TIME, and does not necessarily entail the same 
corresponding reduction in throughput time, since the latter 
is dependent also on processing time, overlapped operations, 
and many other related system parameters. Nevertheless, 
this small example illustrates the power and simplicity of 
this optimization technique. Furthermore, it can be expected 
that the larger the subsystem and the more disparity there is 
among file sizes, record sizes, program run frequencies, and 
memory constraints, the more substantial one can expect the 


gains to be from such a system approach. 
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NONLINEAR PROGRAM TO SOLVE FGR OPTIMAL BLOCKING FACTORS 
FOR TAPE FILES IN A SU8SYSTEM OF TWO OR MORE 
INTERRELATED PROGRAMS. 


FORMAT FOR INPUT DATA: 

CARD 1 (RIGHT-JUSTIFIED INTEGERS) 
COLS Sie J NOe UP IAPE* Five se INSSUBSVSTEM 
COUS 7 6-10 NO- OF CONSTRAINTS 
COLSGL INS LIMIME NOS) OF TTERATIONS 
COES ella 0 START PRINTING FROM THIS ITERATION 
COS SZ 1225 STOP. PRANGING AT THES ITERATION 
EULESS, 26-20 PRINT LTeRATIONSSUP SI HiSSt UT TPLe 


CARDS 2 TO I TNDV LAL PEASI@ UE SOlui ai 

PORG@GC SO PENSE LELDy (Sore LD Sy Pere CAL Oe (Decimal) 
CARDS I+1 TO J NO. OF ALTERNATE TAPE SUFFERS 

2a COOLS? PERC EIELD Ss SUPHELDSE VER? CARD CINTEGER) 
CARDS SAilgirork RECURDESIZES 

LOOCOUS. PERe FIELDS VSL PLELDSE ERs CARDS DECI MAINE 
CARDS K+l TOL PIES SLZES 

LOPCOLSSPERPRIELDs 8 FIELDS PER GARD (DECIMAL)© 
CARDS L+1 TO B FREQUENCIES OF PROGRAM RUNS 

ROZCGEL St PERSIE] Ds Se EITELDS? PEO CARDILUECIMAL) 
CARDS M+1i TON MEMORY CONSTRAINTS MATRIX (BY ROW?) 


ROS COBRG PER FIELD, G FIELDS Pex CALUMUDECINAL) 
IMPLICIT: REAL*38 (A-H,O-Z) - 
DIMENSIGN CL(649)_, VL(9), Y¥C9) - 
DIMENSION FREQ(5)», FSIZE(O}, RSIZE(S), TALT(46) - 


INPUT MATRIX OIMENSTIONS, LIMIT NO. OF ITERATIONS; 
PRINTING PARAHETERS AND CONVERGENCE TOLERANCE. 


READ (5,100) NVyNCoLIMIT, 4FROMsNTOs MUL, EPS 
WRITE (6,L01) NVy NCe LIMIT, EPS 


INPUT INITIAL FEASIBLE SOLUTION 
READ (5,110) (Y¥CI)sT=1lsNV) - | 
INPUT NO. OF ALTERNATE TAPE BUFFERS 
READ (5,710) (IALT(I),I=1,NV) - 


INPUT RECORD SIZES 
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aT ‘BINT a eo oTreesTs ; 
moITIOR ZaGreaws 2e0F 7 
(iaMEINA) GAAS WOT SD IaIT)S eCIN0A anheaeD Oe 
SAD eae SAAT | OCAMART IA AG i ar i: 
(AINOHE) GRAD A949 C4sa0s 4 Jh2399 925" 
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50 


60 


10 
100 
101 


110 
111 
1l2 
201 
204 


205 


206 


READ (59110) -(RSIZE(1),1=1,NV) 
INPUT FILE SIZES 

READ (59110) (FSIZE(I),I=1,NV) 
INPUT RUN FREQUENCIES 

READ (5,110) (FREQ(T),I1=1,NV) - 
INPUT MEMORY CONSTRAINTS MATRIX 


WRITE (6,111) 

DO 50 J=1,NC 

READ (54110) -€((CL(UI_9J),T=1sNV)sVL(J)) 

WRITE S065 112 JUG UEe 3) l=1s NV) Viet) 

WRITE (6,201) (I1yY(I),1=1,NyV) - 

CALL GPROS (ClLrVLe¥ysFaNVsNCyEPSysLIMIT,MFROM,yMTOyMULZIER, 
* FREQ,FSIZE,IT) 

IFCIER.EQ.0) GO TO 60 

WRITE(6,204) ITs(1I,Y(1I),IT=1,NV) 

WRITE UGS 20S) SC Uy TALI (It, LR STZElD sy Ls ESIZECUL) aS 
* IsFREQ(I) s1l=1e,NV) 

WRITE (6,205) F 

Whit Ee(G4250) 71TER 

GO TO 70 

WRITE(6,209) ITslIsV¥(t),T=1lyNV} - 

WRITE(6, 205) (1,ITALT(i)51,RSIZE( 1) gITsFSIZE(I), 
* . IyFREQ(I) y1=13NV) 

WRITE SIG 220 )5F 

FORMAT (4012) 

FORMAT (615,F10.6) 


FGREBAT SE 12475 .OND. CE STAPEVEECES =',14, 
*/* NGO. OF MEMORY CONSTRAINTS ="¢14y 
*/* LIMIT NO. OF ITERATIONS =",14, 
*/* CONVERGENCE TOLERANCE =",E8.1) - 


FORMAT (8F10.2) 

FORMAT (*OMENORY CONSTRAINTS MATRIX (LHS.GEeRHS) 3'9/) 

FORMAT (1X,16F3.0) 

EORMAT (2 01NJ TIAL FEASISLE SSGLUTIONE 
*//,{' CES als = shee) 

FORMAT{(* ONO CONVERGENCE AFTER *,15,° ITERATIONS 3s", 
ef /,(' pg EAS rE) -=",F6.2}) 

FORMAT (ILOAL be wBUEFERS Ste Ll i GRECURDSSLZESs is 
*132,'FILE SIZES:*,150,'RUN FREQs*%, 
xf/f(t ACtei2s*) = 8 Leis eis ce ee) =*,F5.0,132, 
S812, !) =the Us 10s FREQ Us 125 De abies) 

FORMAT('OVALUE OF OBJECTIVE FUNCTION Z =",F9.0) 
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209 FORMAT (*OCONVERGED TO OPTIMAL SCLUTION AFTER *%eI5;, 
ee TVERATIONS Sty / sy UX sce aero ec 


220 FORMAT CTOMINIMUM VALUE OF OBJECTIVE FUNCTION Z =", 


*F9.O) 
250 FORMAT (C*OERROR CODE =* ,J4) - 
TO CONTINUE 

STOP 


END 


SUBROUTINE GPROJ 


PURPOS 


e3 


TO MINIMIZE NONLINEAR OBJECTIVE FUNCTIONS WITH LINEAR 
CONSTRAINTS USING THE GRADIENT PROJECTION METHOD. 


USAGE: 


CALL GPROJ (ChLeVbLy¥sFsNVeNCyEPSyLIMIT,MFROMsMTOsMULy IER, 


FREQ,FSIZE, IT) 


DESCRIPTION OF PARAMETERS: 


CL 


VL 


uf 


LIM 


MFROM 


MTO 
MUL 
IER 


IT 


j 


tof { 


INPUT COEFFICIENT MATRIX (DIMENSION NV*NC) GF THE 

LINEAR CONSTRAINTS IN THE FORM CL*Y.GE.VL ({REAL*8). 

INPUT VECTOR GF R.HeS. OF LINEAR CGNSTRAINTS OF 

DIMNENSICN NC (REAL*8). 

VECTOR OF DIMENSION KV CONTAINS THE INITIAL ARGUMEN 

WHERE THE ALGO21 THE “STR PSS9GN "RETURNS AY GON PAINS 

THE ARGUNENT CORRESPONDING TQ THE COMPUTED WHINIAUM 

FUNCTIGN VALUE {REAL*8). 

CONTAINS MINIMUM FUNCTION VALUE ON RETURN 

(Pete F-Pt yJoe(ReAt*s). 

NUMBER OF VARIABLES CINTEGER*4) « 

NUMBER CF CONSTRAINTS CINTEGER*4). 

TEST VALUE FOR ZERO APPRUXIMATED BY RGUNDOFF ERROR. 

SUDTABEES VALUES “LN VRANGE 9LOEs—4" tO = 104" —5 BUREN 20) « 

MAXIMUM NUMBER OF ITERATIONS. SUGGESTED LIMIT IS 

LIHIT = NC CINTEGER*4). 

START PRINTING FROM THIS ITERATION. 

SICP PRINTING AY Thi Seri eekaArroN. 

PRN Tee ALON Soe 4 aM eter kote 

ERROR CODE (CINTEGER*4). 

TER= 0 ; NO ERROR 

IER= J 3 CONSTRAINT J INVALID OR MISSING 

TER=—T" =" FEASIBLE REGLON UNBOUNDED CICE. NC IS NOT 
GREATER THEN NV) - 

TER=-2 s INITIAL VALUE OF Y oIS NOT FEASISCEE 

TER=-3 =: CONSTRAINTS LINEARLY DEPENDENT (I.E. .Q IS 
GREATER THAN NV) 

IER=-4 3: NC CONVERGENCE IN LIMIT ITERATIONS 

FREQUENCY “EACH TAPE FIEe iS *PRUCESS EUs TNeRUN CYCLE 

OF SUBSYSTeM. (CAN BE FRACTIGNAL) 

AVERAGE NO. CF RECORDS PER -TAPE FILE. 

ITERATIGN CGUNT 

LARGEST VALUE ZUFANVAT OG SSEeM PROCESSED 

NVMAX+1 
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NONDUMMY DIMENSIONS YY( MAX), DF (MAX) »,Z(MAX) yVOQ(MAX) pHEMAX) » 
CQ( MAX¥HAX) »CQT (MAX*MAX) TC MAX#MAX Dy 
TIC MAX*MAX) » RO MAX*MAX ) y PQ( MAX¥MAX) 


ROUTINES AND FUNCTION SUBPROGRAMS REQUIRED: 


MEQAL — DSOS:DSOLIB —- MATRIX EQUATE 

MADD = DS05:DSOLIB? f=" MATRIX? ADD 

MMPY 20505-0500 158 9 —) MAIRIX HULTIPLY (NONDUMMY DIMENSTON 
MINVR -— OSOS:DSOLIB - MATRIX INVERSION (NONDUKMY DIMENSION 
DOT ee oO os USO Tee VCOTOR DU) PRODUCT 

MinAoe — oUPPLIED —- MATRIX TRANSPOSE 

FUNCT -— USER-WRITTEN -— A SUBROUTINE CONCERNING THE FUNCTION 


TOVBES MINT RIZED, MUST BEVOR THE TEORM 
SUBROUTINE FUNCT (€(NsARGsVAL,GRAD) 
AND HUST SERVE THE FOLLOWING PURPOSE: 
TO EACH N-DIMENSIGNAL ARGUMENT VECTGR 
A FUNCTIGN VALUE AND GRADIENT VECTOR 
MUST BE COMPUTED, AND ON RETURN 
SIORED@ING VAL CAND TGRAD  SESPEC IY ELLY 
(NeBe ARGy VAL, GRAD ARE REAL*3). 


REMARKS : 


(1) INPUT MATRIX CL ASSUMED TO BE STORED COLUMNWISE IN 
NYSNE*SUCCES! VE COCATIONS. 

(ZincUTH CLARO V IS ARE VMODIEILED AS (THE. COEFELC LENT SPARE 
NGRMALIZED- IN THE ALGCRITHM. 

(3) NC MUST BE “GREATER THAN NV. 

C4 ENG TTA VALUE GE Y JS REQUIRED SID LE TINT IHe 
PEASIELCS REGION DEF INEDUS Y Thee ULNEARSCUNSTERAINIS .< 

5 ) CAL PSE CUTINESP USED I NY CONUUNCG TIONPHITHM AIS = SUSeRUTINE 
MUST HAVE REAL VARIABLES IN DOUBLE PRECISION FORM 
(Cie. wheal Si. IMIS) CAN CASILY (BE DUNES OYSUS ESO 
THE FORTRAN STATEMENT 

TMPUTC? PY REALS ECA—-H, 0-7) 

AT THE BEGINNING OF EACH ROUTINE} 


REFERENCES: 


(1) KIRK, DeEo,; *CPTIMAL CONTROL THEORY —- AN INTRODUCTION’, 
PRENTICE-HALL, 1970, PP. 373-393. 

(2) ROSEN; J.8-, "THE GRADIENT PROJECTION METHOD FOR 
NONLINEAR PROGRAMMING PART ly LINEAR CONSTRAINTS'» 
JeSlAM, VULe8s 1960s) Pees bi —2Li 
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SUBRGUTINE GPROJ (CL,VL».Y,FsNVsNC-EPS,LINIT, 
MFROM,MTCsMULyITERs FREQ,FSIZE;1IT) 

IMPLICIT REAL*8(A-H,0Q-Z) 

DIMENSION CLIUNVsNC),VLONC) +YCNV) 

DIMENSION FREQUNV), FSIZE(NV) 

DIMENSION YY(7) eDF (7) ¥Z(7) 5 VOUT) sHI7) 9£€Q149) ,CQT(I49), 
T(49),T1(49)9R(49),PQ(49) 

INTEGER¥4 Q 

INTEGER*2 H 

REAL*8 INF/7.D075/ 


INITIALIZE VARIABLES 


IF({NC.LE-NV)GO TO 91 
IT=0 
Q=0 


NORMALIZE LINEAR CONSTRAINTS 


DO 2 J=1;NC 

S=DOTICL (19J),CL0 1,5) _NV) 
LEADABSA SP. LASER SHEEOACL IT 
PE(DAGS (S=1e7 Ul ees) 500102 
S$=1./DSQRT(S) 

DO 1 T=laNV 
CL(IsJJ=S¥*CLUI gd) - 
VL(JI=S*VLOJ) 

CONT INUE 


DETERMINE ACTIVE CONSTRAINTS AND FEASIBILITY CF POINT Y 


DO 5 J=lyNC 

S=DOT(CL (12J)s¥sNV)I-YLEJ) 
IF(DABS(S).-GreEPS)IGO TO 4 
Q=Q+1 

H{Q)=J 

GO TO 5 

IE($-17-0.)G0 TO 92 

CONT INUVE 


DETERMINE PROJECTION MATRIX PQ 
IF(Q.GTeNV+1)G0 TO 93 


IT=17+1 
IFC IT.GT~LIMIT)GO 10 94 


>. 4 Lone #1 1) io 
SMRaTIO. 
an? ty, ae 
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DO 13 I=1,NV 

HG) hs hbeTY 

IJ=14+(J-1)*NV 

Bei GcaP Aa ree) agep she 
PQ(IJ)=0. 

GO TO 13 

PQ(TJ)=1. 

CONT INUE 

IF{Q.£0.0)G0 TQ 15 

DO 14 I1=1,Q 

J=H(T) - 

IJ=14+(I-1)4NV 
CALLSMEOAL(COULJI5CL (14 J) NVe ls 1) 
VQ(I)=VL(J) 

GAELEMIRAS(COT CONV, QO) - 
CALL MMPY(TsCQTsCQ1QyNV_Q) - 
CALL MINVRITI,T DET, Q) 

CALL MHMPY(R,TI,COT29,0,NV) 
CALL MMPY(TsCQsR,NV,Q,NV) 
CALL MADD(PQ,PQsT,NVeNV,—-1) 


CALCULATE GRADIENT PROJECTION Z 


CALL FUNCT(NV,Y,F DF ,FREQsFSIZE) 
CALL MEQAL (DF sDF yNVyly—-1) - 

CALL MMPY(25?Q,0FsNVeNV_1) | 
ZNORM=DSORT(COT(ZsZ_NV)) | 


CALCULATE MATRIX R AND RQ 


RQ=- INF 

IF (@.EQ.9)GO' TO'°18 

CALL MMPY(T+COT,DFsQaNVo Lh} 
CALL MMPY(R.TI19Ts97Q21) 

DO 17 J=1, 
IF({RO-GEe-R(I))GO TO 17 
RQ=R{(1) 

K=] 

CONTINUE 


TEST FOR CONSTRAINED MINIMUM. 


IF(ZNORM.GT.EPS)IGO TO 20 
IF(RQ.LE.EPS)GO TO 98 
GO TO 25 


CALCULATE BETA 
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IF(Q.EQ.0)GO TO 26 
BETA=0. 

DOM227= 150 

S=0. 

DO 21 J=I,Q 
IJ=14+(J-1)*Q 
S=otVADSt Lit ral) ) 

PER iAe Liss) CeLA=S 
CONTINUE 
IF(RQ.LE-BETA)GO TO 26 


DROP HYPERPLANE H(K) FROM ACTIVE CONSTRAINTS 


H(K)=H(Q) 
Q=Q-1 
GO TO ill 


FORM NORMALIZED GRADIENT PROJECTION Z 


S=1e /ZNORM 
DO 27 I=1;NV 
ZCI) =S*Z01) » 


DETERMINE MAXIMUM POSSIBLE STEP SIZE (STEPM) 


STEPM=INF 

00 30 J=1,NC 

IF (Q.6Q@.0)G60° TO 29 

DO 28 T=1;Q 

Petje Oen CL) 6010 o30 

CONTINUE 

S=DOT{ICL (lyJ) sZ7NV) 
LF(DABS(S).-LT-eEPS)GO TO 30 
SFEC=UVL CI) —O06T (CELTS SIs YeNVIIIS 
TE(STEPS LE. UeCRSSTEP<GiesTEPMIGO 107530 
STEPM=STEP 

M=J 

CONT INUE 


DETERMINE TENTATIVE NEXT POINT YY 


DO 36 IT=1ysNV 

YY (1T)=Y(1)4+STEPM*Z(1) | 

CALL FUNCT(NV,YY¥eF,DFsFREQsFSIZE) 
CALL MEQAL(DFsDF yNV_1l_-1) 

IE (DOT(DFE »ZsNV).GE.-EPS)GO TO 45 


FIND NEXT POINT YY BY BISECTION SEARCH 


80 


2 


i] 


40 TH=STEPM/2. 
STEPM=TH 
41 DO 42 I=1;NV 
42 YY(I)=Y (I) +STEPM*Z(1) 
CALL FUNCTINV,YY,F;0F,FREQ,FSIZE) 
CALL MEQAL{DF»DFyNV,1,-1) 
S=DOT (OF aiyNV) 
IF (DABS(S).LT~EPS)GO TO 46 


TH=TH/2. 
STEPH=STEPM+OSIGN(1.00,S)*TH 
GO TO 41 
C 
C ADD HYPERPLANE M TO ACTIVE CGNSTRAINTS 
C 
45° Q=Q+1 
H(Q) = 
C 
C UPDATE Y AND START NEXT ITERATION 
C 
46 CALL MEQAL(YsYY¥eNVo1ly91) - 
c . 
C PRINT INTERMEDIATE ITERATIONS 
C 


ML=IT/MUL*AUL 
TFECIT «LT eMEROM. OR IT. GPeHTO.OR.IT.NEsRL) -GO 10°10 
RRITE (6,1000) IT, (19 YC(I)sI=1,NV) 
1000 FGRHAT)'( *OSOLUTION AFYER ~*,15,° PTERAT IONS ss bs 
xf fi" Keel =e Gor) ) 


GO TO 10 
C 
C RETURN ERROR CODE 
C 
91 IER=-1 
GO TO 99 
92 IER=~2 
GO TO 99 
93 IER=-3 
GO TO 99 
94 IER=-4 
GO TO 99 
97 IER=J 
GO TO 99 
98 IER=0 
99 RETURN 


END 
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SUBROUTINE MADD (AsByCsNRyNCeKS) 


IF KS.LT.0, SUBROUTINE MADD SETS A EQUAL TO B. 
IF KSeGE.0; SUBROUTINE MADD SETS A EQUAL TO B PLUS Ce 


IMPLICIT REAL*38(A-H,0-Z) 

DIMENSION ACNR NC) s8(0NRgNC) sC(NRNC) 
IF(KS.GE.0) GG TO 200 

DO 150 I=1,NR 

DO 150 J=1,NC 

AC(IsJ)=B(1,J)I-ClI,J) 

RETURN 

DO 250 I=1,NR 

DO 250 J=1;NC 


250 A(I,J)=68(1;J3)4+C(1I2J) 


RETURN 
END 


Ck ee 


150 


200 


250 
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SUBROUTINE MEQAL (A;B,NR,NC,KS) 


IF KSeLT-Os SUBROUTINE MEQAL SETS A EQUAL TO MINUS B. 
IF KS.GEeO, SUBROUTINE MEQAL SETS A EQUAL TO Be 


IMPLICIT REAL*8( A-H,0-Z) 
DIMENSION A(NR NC) ,B(NRSNC) 
IF(KS.GE.0) GO TO 200 
DO 150 I=1,NR 

DO 150 J=1,NC 
AtI,J)=-Bl(I2J) 

RETURN 

DO 250 FT=l1;sNR 

DO 250 J=1;:NC 
AlCI,J)=B(1I,J) - 

RETURN 

END 
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SUBROUTINE MMPY (A,ByCyNRAyNCB;NCA) 


SUBROUTINE MMPY SETS MATRIX A EQUAL TO MATRIX B 
POS iarnUe ee PCO Yea TALIA cure 


NRA = NUMBER OF ROWS IN AyB 
NCB = NUMBER OF COLUMNS IN 8, ROWS IN C 
NCA = NUMBER OF COLUMNS IN A,C 


IMPLICIT REAL*8(A-H,O-Z) 
DIMENSIGN A(1024} ,811024) ,C(1024) ,D( 1024) - 
DO 100 I=l1~sNRA 

DO 100 J=lyNCA 

MD={(J-1) *NRA+I 
D{MD)=0.0 

DO 100 K=1,NCB 

MA=(K-1) *NRA+I 
MB=(J-1)*NCB+K 
D(MD)=D(MO)4+B( MA) *¥COMB) 
NME=NRA*NCA 

DO 200 JT=1;NMNE 

A(1I)=D{tI1) 

RETURN 

END 
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SUBROUTINE MTRAS (RyAyMyN) 


SUBROUTINE NTRAS TRANSPOSES MATRIX Ae 


IMPLICIT SREAL4S(A—H, 0-2) 
DIMENSION ACM sN) »ROUN,M) 
DO 1 I=lsM 

DO 1 J=1:N 

RiJdsI)=A(TsJ) 

RETURN 

END 
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SUBROUTINE MINVR (AyBsDET,N) 


SUBROUTINE MINVR COMPUTES THE INVERSE AND DETERMINANT 
OF MATREX Bs 


IMPLICIT REAL*8{ A-H,O0-Z) 
DIMENSION A(NeN) »B(NgN) 
DIMENSION PIVOT(32),I1PIVOT(32),INDEX(32,2) 
EQUIVALENCE (IROWsJROW) s CICOLUM, JCOLUM) » CAMAXa Te SHAP) - 


INITIALIZE VARIABLES 


DET=1.0 

00 40 I=1y,N 
IPIVOT{1I)=0 
DO 40 J=l1sN 
A(I,J)=B(1eJ) 


PERFORM INVERSION ROUTINE 
DO 600 I=1;5N 
SEARCH FOR PIVOT ELEMENT PIVOT(I)« 


AMAX=0.0 

DO 110 J=1;N 

IFCIPIVOT(J).EQ9.1) GO TO 110 

DO 100 K=15N 

IF (CIPIVOT(K)—-1) 80;100;1000 

IF (DA3S(ANAX) -GEDABSLATJ,K))) GO TO 100 
TROW=J 

ICGLUM=K 

AMAX=A(J 9K) 

CONTINUE 

CONTINUE 

IPIVOT (ICGLUM)=1P1TVGT(ICOLUM)+1 


INTERCHANGE ROWS TO PLACE PIVOT ELEMENTS ON DIAGONAL. 


IF(CIROW.~EQ-ICCLUM) GO TO 260 
DET=—-DET 

DO 200 L=1;N 

SWAP=A(ITROW,L) 
A(IROW,L)=ACICOLUMeL) 


200 AC{ICOLUM,L)=SWAP 
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260 


350 


450 
550 
600 


T00 
710 
LOCO 


INDEX(1,1)=I1ROW 
INDEX(1T,2)=I COLUM 
PIVOT(1I)=ACICOLUM, ICOLUM) 
DET=DET*PIVOT(I) 


NORMALIZE PIVOT RGW CN PIVOT ELEMENT. 


A(ICOLUH ,ICGLUM)=1.0 
DO 350 L=1,N 
ACICOLUM sL J=ACICOLUM,L)/PIVOT(I) 


REDUCE NON-PIVOT ROWS. 


DO 550 L1=1,N 
IF(L1.EQ.1CGLUM) GO TO 550 
T=A(L1,1COLUN) 
A(L1,I1COLUM)=0.0 

DO 450 L=1;N 
A(L1sL)=A(0L1,L)-ACICOLUM,L)¥T 
CONTINUE 

CONTINUE 


INTERCHANGE COLUMNS. 


DO 710 I=1,N 

L=N+1-I 

IF (INDEX(Ls1).EQS-INDEX(L22)) GO TO 710 
JROW=INDEX(Ls1) 
JCOLUM=INDEX(L,2) 

DG 700 K=i;N 
SWAP=A(K,JROW) 

AK, JROW) =ACK,JICOLUM) 
AK, JCOLUM} =SWwAP 
CONTINUE 

CONTINUE 

RETURN 

END 
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FUNCTION DOT (X+YsNN) 


SUNCT IONS DOTS CALCULATES THE DOUI—PRODUCTS OF 
NN-DIMENSICNAL VECTORS X AND Ye 


TIPE IGT SREAU-es CA-H, 0-2) 
DIMENSION X(NN) sYCNN) 
DPR=0.0 

DO 100 I=1,NN 
DPR=DPR+X(1)*Y(T) 
DOT=DPR 

RETURN 

END 
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SUBROUTINE FUNCT (NyXyVALsGRAD, FREQ, FSIZE) - 


SUBROUTINE FUNCT COMPUTES THE VALUE GF THE 
OBJECTIVE, FUNCTION ANDaLTS PARTIAL DERIVATIVES 
WITH RESPECT 10 EACH VARIABLE AT POINT X- 


IMPLICIT REAL*8 (A-H, 0-Z) 

DIMENSIGN X(N),5y GRAD(N), FREQ(IN), FSIZE(N) 
VAL=0- 

DO 10 I=1;N 
VAL=VAL+FREQ(I)*FSIZECIT) ZX(1) - 

GRAD (1 )=-FREQ(1)*FSIZE(1) /XC1) ¥¥*2 

RETURN 

END 
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(LHS.GEeRHS) 3 


—-140. O. 
0. ~2700. 
0. coe ORB 
Oe 0. 
Oe 0. 
0. O. 
le 06. 
0. le 
0. 0. 


SAMPLE OUTPUT : 

NOs Ut a APESEILES 

NO. OF MEMORY CONSTRAINTS 

LIMIT NO. OF ITERATIGNS 

CONVERGENCE TOLERANCE 

MEMORY CONSTRAINTS MATRIX 

~300. —-300. = TOs 
~300. 0. 0. 

O. Oe 0. 
le 0< 0. 
Ce. le 0. 
0. O. le 
0. Oe O« 
Oe 0. 0. 
O. 0. O. 

INITIAL FEASIBLE SOLUTION 

a ie ieee Geral) 

X€ 2) = 10.00 

X( 3) = 10.00 

X€ 4) = 10.00 

X{ 5) = 10.00 

X( 6) = 10.00 

SCLUTION AFTER 100 ITERATIONS : 

teal —iciee 

X(6 20 Si 2e25 

Xe sae wl 

X€ 4) = 10.89 

X{( 5) = 24.43 

X( 6) = 10.76 

SOLUTION AFTER 200 ITERATIONS =: 

X¢( 1) = 12.24 

X(72) “= 12524 

X{ 3) = 7.214 

X{ 4) = 10.80 

X¢ 5) = 34.272 

X{ 6) = 54228 

SOLUTION AFTER 300 ITERATICNS =: 

X{ 1) = 12.24 

X( 2) = 12.24 

at Sep ee 

X( 4) = 10.30 

X( 5) = 25.03 

X( 6) = 4297 


‘2 @alranseD 604 
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SOLUTION AFTER 400 ITERATIONS =: 


X( 1) = 12.24 
Wee esol a ot 
X( 3) = 7.214 
X€ 4) = 10.80 


OR Ba 35216 


NO CONVERGENCE AFTER 501 ITERATIONS ; 


X€e1) = 12.24 

X( 2) = 12.24 

Kes) = 11s 

X( 4) = 10.80 

AAS) = 635621 

X{ 6) = 4019 

Aisle OU ELK < RC OURD) 512 > ceils bres. RUN FREQ: 
At 1) = 2 RC P)e= 150. Piel) 2-2 0000. FREQ( 1) 
At 2) = 2 R{ 2) = 150. El 2)i= 350000. EREOL 2) 
A( 3) = 2 R{ 3) = 80. F{ 3) = 105000. PREQtUS) 
AC 4) = 2 R{ 4) = 70-6 F{ 4) = 20000. FREQ( %) 
(REC SY fe 74 Rif 5b Je 1006 F( 5) = 5000. PREOQ(SIS) 
RIS6) = 2 R¢ 6) = 100. F{ 6) = 500. FREQ{ 6) 


VALUE CF OBJECTIVE FUNCTION Z= 270516. 


ERROR CODEE= W=4 
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